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14.  ABSTRACT 


The  atmospheric  transmission  windows  of  2-5  and  8-12  pm,  coupled  with  organic  and  other  chemical  absoiption 
lines  occurring  throughout  this  middle-infrared  (mid-IR)  wavelength  region  give  rise  to  a  wide  variety  of  medical, 
scientific,  commercial  and  military  applications.  Communications,  remote  sensing,  IR  countermeasures,  laser 
surgery  and  non-invasive  imaging  are  just  a  few  of  the  drivers  of  high-power  solid-state  mid-IR  laser  development. 
These  laser  sources  must  be  versatile  enough  to  operate  in  a  variety  of  temporal  modes  from  continuous  wave  (CW) 
all  the  way  to  ultrashort  pulse  while  still  being  widely  tunable  for  wavelength  agility.  All  of  this  is  required  at  ever 
increasing  power  output  levels  while  conforming  to  size,  weight  and  power  consumption  limitations  under  harsh 
operating  environmental  conditions. 

Chromium-doped  zinc  selenide  (Cr2+:ZnSe)  lasers  operating  in  the  2-3  pm  region  are  excellent  candidates  to  help  fill 
these  vital  roles.  As  a  transition-metal  doped  II- VI  chalcogenide,  Cr2+:ZnSe  has  a  number  of  positive  advantages 
over  existing  laser  sources.  Development  and  power  scaling  of  these  lasers  however,  has  been  hampered  by  thermal 
issues  which  have  so  far  limited  the  ability  of  these  lasers  to  be  applied  to  systems-level  development. 

This  work  presents  research  into  the  nature  and  mitigation  of  these  critical  thermal  issues  in  development  of  versatile 
Cr2+:ZnSe  laser  sources.  Advanced  models  for  thermal  and  laser  performance  are  developed  and  used  to  design 
optimally  configured  laser  systems.  Among  other  advances  for  this  material,  >10  W  CW  output  from  a  Cr2+:ZnSe 
oscillator  and  master-oscillator  /  power  amplifier  systems  producing  multi-watt,  widely  tunable  power  levels  are 
demonstrated. 

Advisor:  Peter  E.  Powers,  PhD 
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Chapter  1 
Introduction 


1.1  Aims  of  this  Work 

Laser  sources  operating  in  the  middle  infrared  (mid-IR)  wavelength  region  from  2-12  pm  are 
necessary  for  a  number  of  military,  scientific  and  other  types  of  applications.  These  applications 
are  driven  in  general  by  both  the  atmospheric  transmission  bands  of  2-5  pm  and  8-12  pm  as  well 
as  the  locations  of  chemical  lines  throughout  the  region.  The  very  nature  of  applications  such  as 
remote  chemical  sensing,  infrared  countermeasures  (IRCM),  laser  communications  and  laser 
surgery  drive  other  laser  operational  requirements  such  as  output  power,  tunability,  beam  quality 
and  room  temperature  operations  to  name  just  a  few.  The  broad  goal  of  this  research  was  to 
understand,  design  and  construct  useful  laser  sources  with  emission  in  this  important  wavelength 
range.  This  work  builds  on  years  of  basic  research  in  this  area  but  also  incorporates 
requirements-driven  feedback  from  systems-level  developments  with  the  end  goal  of 
accomplishing  a  broader  understanding  of  the  science  involved  while  still  producing  useful  end 
products. 

The  approach  used  to  produce  coverage  of  these  broad  bands  was  to  develop  high  power,  tunable 
laser  direct  laser  sources  to  cover  the  2-3  pm  region  and  then  use  nonlinear  devices  to  convert 
this  power  out  into  3-5  pm  and  8-12  pm.  This  work  focuses  on  the  first  piece  of  that  puzzle: 
developing  versatile  direct  laser  sources  in  the  2-3  pm  wavelength  range.  Modeling  of  system 
behavior  was  used  whenever  possible  to  save  research  time  and  enhance  the  probability  of 
success  by  prioritizing  research  alternatives  based  on  model  results. 

1.2  Dissertation  Outline  and  Relevant  Publications  and  Presentations 

This  dissertation  is  broadly  broken  up  into  chapters  covering  different  types  of  CL  :ZnSe  lasers 
which  were  investigated  under  the  umbrella  of  this  work.  Chapter  2  will  cover  the  background 
of  high-power  mid-IR  lasers  as  well  as  a  review  of  relevant  research  and  our  previous  work  in 
this  area.  Chapter  3  presents  modeling  of  thin-disk  lasers  undertaken  in  an  attempt  to  better 
understand  the  effect  of  different  pumping  sources,  both  CW  and  pulsed,  on  the  thermal  load  of 
the  material.  Chapter  4  covers  work  done  to  develop  a  Cr2':ZnSe  fiber  laser.  Chapter  5 
showcases  work  done  to  model  and  build  a  high  power,  tunable,  thin-slab  Cr  :ZnSe  master 
oscillator  /  power  amplifier  system.  The  publications  and  presentation  associated  with  each 
chapter  can  be  seen  in  Table  1. 


Table  1:  Dissertation  related  presentations  and  publications 


Chapter 

Presentations  and/or  Publications 

1  Introduction 

2  Background  and  Motivation 

SPIE  Newsroom  Article  [l]b 

3  Time-dependent  Modeling 

IEEE  Journal  of  Selected  Topics  in  Quantum 
Electronics  [2]b,  2006  Phot.  West  presentation3, 
SPIE  Proceedings  [3]a 

4  Cr  .ZnSe  Fiber  Laser 

5  CL  :ZnSc  High-power  Systems 

a  t-»  .  i  .  i  b  .  i 

2009  ASSP  presentation,  [4],  2010  Phot.  West 
presentation,  SPIE  proceedings  [5]a,  2010  Optics 
Express  [6]a’°,  2010  CLEO  presentation  [7]a,c 

Cl  1 

a  Presenter  and/or  primary  author; b  co-author; c  submitted 
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Chapter  2 

Background  and  Motivation 


2.1  Introduction 

As  mentioned  in  Chapter  1,  the  2-12  pm  wavelength  range  is  one  in  which  a  number  of 
important  applications  occur.  Due  to  the  absorption  characteristics  of  its  constituent  gasses,  the 
atmosphere  has  two  broad  ranges  from  2-5  pm  and  8-12  pm  in  which  high  transmission  occurs 
(Figure  1).  Applications  such  as  laser  ranging  and  imaging,  target  identification  and  laser-based 
communications  which  rely  upon  sending  and  receiving  laser  pulses  leverage  these  high- 
transmission  windows  in  order  to  extend  operating  distances.  Anti-aircraft  missiles  utilize  these 
windows  to  maximize  their  ability  to  detect  heat  from  target  aircraft  either  generated  by  engine 
combustion  or  by  friction  between  the  aircraft  and  the  atmosphere.  This  threat  drives  the 
requirements  for  lasers  in  these  ranges  for  use  in  infrared  countermeasures  for  aircraft  self¬ 
protection. 


Wavlength  (microns) 


Figure  1:  Simulated  atmospheric  transmission  indicating  areas  of  interest  for  mid-IR 
laser  development 

In  addition  to  the  atmospheric  absorption  characteristics,  there  are  a  number  of  organic  and  other 
molecular  absorptions  which  occur  in  this  region.  These  absorption  features  drive  applications 
such  as  remote  chemical  sensing  in  which  lasers  are  used  to  probe  areas  from  a  safe  or  practical 
distance  to  determine  the  makeup  of  unknown  samples.  Similar  absorption  properties  of  tissues 
and  proteins  drives  medical  applications  such  as  laser  tissue  cutting,  ablation  and  welding  [8,  9]. 

All  of  these  applications  have  specific  laser  operating  requirements.  Temporally,  the  laser  may 
need  to  operate  anywhere  from  continuous  wave  (CW),  to  long  pulse  (ms),  to  short  pulse  (ns),  all 
the  way  down  to  ultrashort  pulse  (ps,  fs).  Many  of  the  applications  require  the  ability  to  tune  the 
wavelength  through  tens  of  nanometers  or  more  in  short  time  periods,  allow  for  frequency 
modulation  or  multiplexing  or  select  specific  absorption  features.  The  end  result  of  all  of  the 
above  is  a  broad  requirements  field  which  necessitates  multiple  lasers  to  ensure  coverage.  What 
is  really  needed  is  a  laser  versatile  enough  to  do  the  work  of  many  lasers. 

2.2  Early  Solid-state  Lasers 

Townes  and  Schawlow  laid  the  groundwork  for  the  first  solid  state  laser  in  1958  with  their 
description  of  infrared  and  optical  masers  [10],  which  lead  to  Maiman’s  demonstration  of  the 
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first  ruby  (chromium-doped  sapphire)  laser  in  1960  [11].  Attention  quickly  turned  to  rare  earth- 
doped  materials  including  Nd:YAG  [12],  Nd:YALO  [13],  among  many  others.  Now  Nd-based 
lasers  are  ubiquitous  with  output  at  1.064  pm  along  with  many  possible  nonlinear-converted 
wavelengths.  These  same  hosts  have  also  been  demonstrated  with  a  wide  variety  of  other  rare 
earths  such  as  erbium  [14],  and  holmium  [15]  just  to  name  a  few  [16].  Rare  earth  lasers  fill  a 
number  of  vital  roles  in  commercial,  medical  and  military  needs;  however  their  tunability  and 
operating  wavelengths  in  the  mid-IR  range  of  interest  are  limited  by  their  low  crystal  field 
interactions  with  their  host  materials[17].  This  arises  from  the  fact  that  the/ orbitals  in  which  the 
valence  electrons  lie  are  shielded  from  host  interactions  by  filled  outer  orbitals  [18,  19]. 

2.3  Transition  Metal  Lasers 

Transition  metal  (TM)  lasers,  in  contrast  to  the  rare  earths,  have  always  been  found  to  have  broad 
tunability  (with  the  exception  of  the  ruby  laser).  This  comes  from  the  fact  that  the  active 
electrons  in  these  ions  he  in  partially-filled  d  orbitals  which  are  not  shielded  by  filled  outer 
shells.  This  allows  these  ions  to  interact  strongly  with  their  host  crystal  lattices  giving  broad 
emission  and  absorption  spectra.  Early  TM  lasers  were  demonstrated  with  Ni2  [20]  and  Co2+ 
[21]  but  the  tunability  of  these  lasers  was  limited  by  excited  state  absorption  (ESA)  and  were 
only  ran  efficiently  at  cryogenic  temperatures  [22-24],  a  requirement  limiting  the  size,  weight 
and  power  of  the  lasers.  This  inability  to  run  efficiently  at  room  temperatures  is  due  to  the 
decreasing  fluorescence  lifetime  associated  with  increasing  non-radiative  transition  rates  at 
higher  temperatures. 

The  real  breakthrough  in  mid-IR  TM  lasers  came  in  1995  when  a  whole  new  class  of  TM-doped 
II- VI  chalcogenide  materials  was  reported  by  researchers  at  Lawrence  Livermore  National 
Laboratory  [25-27].  Room  temperature  lasing  has  been  reported  in  a  number  of  different 
transition-metal  doped  II-VI  materials  including  Cr2+:ZnSe  [27-29],  Cr2+:ZnS  [27,  28,  30], 
Cr2+:Cdi_xMnxTe  [31,  32],  and  Cr2+:CdSe  [33,  34]  crystals.  Much  of  the  focus  however,  has 
been  on  Cr2  :ZnSe  which  has  the  desirable  properties  of  broad  tunability  [35,  36]  and  absorption 
bands,  a  large  gain  cross  section,  no  excited  state  absorption,  and  finally  high  thennal 
conductivity,  infrared  transparency  and  quantum  efficiency  at  room  temperature.  The  versatility 
of  this  laser  material  has  been  shown  through  fiber  laser  pumping  [37],  diode  pumping  [38], 
continuous  wave  (CW)  lasing  with  efficiency  exceeding  60%  [39,  40],  output  power  levels 
exceeding  12  W  [4],  gain  switched  lasing  with  average  output  power  of  up  to  18.5  W  [41],  as 
well  as  narrow-linewidth  operation  [42].  Only  relatively  recently  have  higher-quality  Cr2  :ZnS 
crystals  been  fabricated  and  because  of  this,  these  types  of  lasers  have  not  been  as  extensively 
demonstrated.  However,  they  have  been  shown  to  be  capable  of  multi-watt  output  [43]  and  can 
have  some  advantages  over  Cr"  :ZnSe  with  respect  to  physical  properties  [44,  45].  Power 
scaling  of  chromium  lasers  has  long  been  hampered  by  the  problem  of  thermal  lensing  due  to  the 
high  thenno-optic  coefficients,  compared  to  other  materials  in  Table  2,  of  ZnS  (46  x  10'6  1/K), 
ZnSe  (70  x  10'6  1/K),  CdSe  and  CdTe  (both  -100  x  10'6  1/K  )  [46].  Thennal  lensing  in  the  gain 
material  leads  to  optical  damage  or  cavity  instability  which  limits  the  maximum  achievable 
output  power.  Thermal  lensing  occurs  when  variations  in  the  refractive  index  across  the 
material,  arising  from  thermal  gradients,  alter  the  phase  profile  of  light  passing  through  the 
material,  effectively  acting  as  a  lens.  This  can  cause  cavity  instability  as  well  as  material  damage 
from  focusing  at  high  powers.  The  thermal  gradients  that  cause  thennal  lensing  typically  come 
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from  excess  heat  generated  by  pumping,  either  through  the  quantum  defect  or  from  non-radiative 
losses  [2]. 


Table  2:  Material  properties  of  various  laser  host  materials [46] 


Property 

ZnSe 

ZnS 

CdSe 

CdTe 

Sapphire  (AI2O3) 

k  (Win'lC1) 

18 

17 

4 

7.5 

27 

dn/dT(10-6/K) 

70 

46 

98 

100 

12 

n  (~A,ias;ng) 

2.45 

2.27 

2.47 

2.7 

1.76 

Hardness 

(Knoop) 

160 

120 

70 

45 

1800,  2200 

These  thermal  issues  can  be  addressed  through  appropriate  choices  of  crystal  geometry,  a  longer 
wavelength  pump  source  and  heat  removal,  all  of  which  help  to  reduce  the  thennal  load  in  the 
gain  element  [2].  A  number  of  designs  for  thermal  lens  mitigation  resonators  are  found  in  the 
literature  and  include  the  fiber  waveguide-based  laser,  ziz-zag  mode  slab  laser  and  the  thin  disk 
laser.  While  some  progress  has  been  made  in  addressing  this  problem  using  thin  disk  laser 
designs,  there  is  still  room  for  improvement  [47]. 
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Chapter  3 

Modeling  of  Time-Dependent  Thermal  Effects 


In  order  to  design  lasers  which  mitigate  the  thennal  effects  arising  from  thermal  lensing  in  ZnSe, 
we  must  understand  the  time-evolution  of  the  lensing.  As  a  first  step  in  this,  a  modeling  system 
was  developed  that  could  be  used  to  gain  an  understanding  of  the  time-dependence  of  thennal 
lensing  as  well  as  to  help  expedite  the  development  of  better  laser  designs.  With  this  in  mind,  it 
was  important  to  not  only  construct  models  with  varying  time-dependent  variables  but  to  also 
study  the  effects  of  variations  in  the  ancillary  model  parameters.  These  extra  steps  were  needed 
to  decouple  the  modeling  assumptions  from  the  parameters  to  be  tested  and  establish  a  baseline 
model  from  which  the  proposed  designs  can  be  constructed. 

When  modeling  heat  distribution  and  material  properties,  finite  element  analysis  is  the  modeling 
method  of  choice.  Commercially  available  finite  element  software  was  used  to  model  the  time 
dependence  of  thermal  lensing  and  temperature  rise  in  a  Cr2+:ZnSe  thin  disk.  Two  general  types 
of  sources  are  available  to  use  as  a  pump  for  Cr  :ZnSe  lasers,  diode-pumped,  Q-switched  lasers 
such  as  Tm:YLF  or  TrmHoYLF,  or  diode-pumped  CW  fiber  lasers  such  as  Timfiber,  which  can 
be  modulated  to  allow  for  varying  duty  cycles.  Both  of  these  cases,  chopped  CW  and  Q- 
switched  pumping,  were  considered. 

3.1  Methodology 

3.1.1  General  assumptions 

To  simplify  the  modeling  process,  certain  assumptions  were  made.  The  laser  source  was 
assumed  to  be  uniformly  absorbed  along  the  depth  of  the  laser  sample.  In  a  typical  single-pass 
absorption,  one  would  expect  pump  absorption  which  followed  Beer’s  law  along  the  depth  axis. 
In  the  case  which  we  have  modeled,  the  pump  is  reflected  at  the  rear  face.  This,  along  with 
multiple  pump  passes  (16  in  our  model)  will  even  out  the  absorption  at  different  sample  depths. 
Heat  transfer  to  the  heat  sink  was  assumed  to  experience  no  losses  due  to  the  interface  and  the 
heat  sinks  maintained  a  constant  temperature  on  their  surface.  Lastly,  conductive  cooling  by  air 
was  neglected  such  that  all  non-heat  sink  surfaces  were  considered  to  be  insulating.  While  there 
will  be  some  conductive  heat  loss  to  the  air,  the  conductivity  of  air  is  orders  of  magnitude  lower 
than  copper  (0.024  WoT’K'1  vs  401  Wm^K"1)^]. 


The  pump  beam  was  assumed  to  have  a  Gaussian  transverse  intensity  beam  profile,  given  by 


fn(x,y)  =  exp 


2((x-x0)2  +  (y  -  y0  )2 )” 


(3.1) 


where  x0  and  y0  define  the  center  point  of  the  beam,  n  is  the  Gaussian  order,  and  w  is  the  beam 
radius  at  which  the  value  of  the  function  drops  to  1/e".  The  Gaussian  order,  n,  could  be  varied  to 
study  the  effects  of  different  beam  profiles,  as  seen  in  fig.  In  the  case  of  our  pump  lasers 
however,  n  =  1  was  used  throughout. 
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An  integration  factor  for  each  of  these  profiles  had  to  be  calculated  to  normalize  each  one  to  1  W 
of  absorbed  power.  Custom  MatLab  functions  (see  Appendix)  were  used  to  calculate  the  heat 
source  (Gaussian)  and  heat  source  integration  constant  for  the  beam  profile.  Two  additional 
custom  functions  were  used  to  verify  that  the  integration  constants  were  essentially  the  same  for 
cylindrical  co-ordinates.  Using  the  results  of  these  calculations,  the  integration  constants  were 
applied  in  the  model  to  normalize  the  power  to  1  W. 


3.1.2  Model  Construction 

The  models  were  constructed  using  FemLab®  3.0a,  version  number  3.0.0.228.  The  first  step  in 
the  model  construction  process  was  to  create  the  model  geometry  using  the  program’s  CAD 
interface.  In  this  step  it  is  important  to  faithfully  recreate  the  important  elements  of  the  design 
while  minimizing  the  extraneous  details.  This  step  is  important  due  to  the  fact  that  the 
complexity  of  the  model  and  therefore  the  amount  of  time  and  memory  required  in  solving  the 
model  are  directly  related  to  the  size  of  the  model  and  the  number  of  surfaces. 

3. 1.2.1  Mode I  design 

A  common  laser  design,  the  face-cooled  thin  disk,  was  simplified  to  consist  of  the  laser  material 
and  a  heat  sink.  Due  to  the  general  assumption  of  perfect  heat  sinks,  the  only  part  of  the  heat 
sink  that  needed  to  be  included  was  the  surface  that  was  in  direct  contact  with  the  laser  material 
or  window.  To  illustrate  this,  Figure  3  shows  how  the  face-cooled  thin  disk  design  was 
‘translated’  into  the  model  geometry. 
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Figure  3:  Basic  face  cooled  design  translation  to  model  geometry 


The  next  step  in  the  model  construction  is  to  define  the  boundary  conditions  for  the  geometry.  In 
general,  there  are  three  types  of  boundaries  in  this  type  of  model.  Any  surface  that  touches  air  is 
considered  to  be  an  insulator  such  that  any  radiative  cooling  by  air  is  neglected.  Any  surface  that 
is  to  be  considered  a  heat  sink  or  that  touches  coolant  should  be  set  to  the  value  of  the  heat  sink 
temperature  constant.  All  other  surfaces  are  usually  intersections  of  objects  or  subdomains  and 
their  parameters  are  solved  for  by  the  program. 

Once  the  model  geometry  was  established,  boundary  conditions  and  subdomain  material 
parameters  were  set  according  to  the  material  properties  of  the  system.  These  parameters  include 
the  important  physical  properties  of  the  host  material,  zinc  selenide.  The  subdomain  settings  also 
include  the  heat  source  Q.  This  is  an  extremely  important  parameter  as  it  not  only  includes  the 
beam  profile  from  above  but  also  the  time  dependence  of  the  pump  beam. 

3. 1.2. 2  Time- depen  dent  pump  algorithm 

A  method  for  simulating  periodic  pumping  by  a  laser  source  was  developed  for  use  in  the  heat 
source.  The  variable  t  is  available  for  use  with  the  time-dependent  solver  and  can  be  used  to 
create  a  periodic  square  wave  pulse.  This  function  had  to  have  the  value  of  1  when  the  pump 
laser  is  on  and  0  when  it  is  off.  The  simplest  method  to  create  a  square  wave  was  to  use  logic 
functions  to  turn  the  power  on  and  off.  However,  solutions  generated  with  this  method  were 
found  to  be  extremely  inconsistent  due  to  discontinuities  of  first  derivates  at  the  on/off  transition 
times.  A  better  way  to  create  a  periodic  square  wave  was  through  the  use  of  a  smoothed  step 
function  known  as  a  smoothed  Heaviside  function.  This  function,  while  still  essentially  the  unit 
step  function  (x  >=  0),  is  smoothed  at  the  step  allowing  it  to  have  a  first  derivative  at  that  point. 
The  Heaviside  function  is  not  a  periodic  function  however,  so  it  needed  to  operate  on  a  periodic 
function  to  achieve  the  desired  pulse  train.  To  understand  the  algorithm  better,  it  is  necessary  to 
look  at  how  a  Heaviside  function  would  behave  when  operating  on  a  periodic  function  such  as  a 
sine  wave  (Figure  4a).  The  Heaviside  function  will  then  turn  on  when  sin  ( 2nt/T)  >=  0  and  off 
when  sin  (2nt/T)  <  0  (Figure  4b)  where  T  is  the  period  of  the  sine  function.  This  produces  a 
square  wave  pulse  with  a  50%  duty  cycle  because  a  sine  function  spends  50%  of  the  time  above 
0  and  50%  below. 
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Figure  4:  Generation  of  a  time  dependent  square  wave,  (a)  sine  function,  (b)  50%  duty  cycle  smoothed 
Heaviside  function,  (c)  sine  function  offset  by  cos  (0.2n)  =  0.628,  (d)  20%  duty  cycle  smoothed  Heaviside 
function. 


To  model  a  square  wave  pulse  with  a  different  duty  cycle,  the  sine  function  that  the  Heaviside 
function  operates  on  must  be  shifted  up  or  down  by  an  offset  so  that  it  spends  the  desired  duty 
cycle  above  zero.  To  determine  the  correct  amount  of  offset,  we  simply  take  the  value  of  a 
cosine  at  the  desired  pulse  width.  For  example,  if  we  want  a  20%  duty  cycle  we  shift  the  sine 
function  by  cos  (0.27t)  (Figure  4c).  This  produces  the  following  square  wave  pulse  with  a  20% 
duty  cycle  (Figure  4d).  Note  that  the  pulse  is  centered  around  the  ti/2  point.  This  algorithm  can 
be  repeated  indefinitely  to  produce  a  time-dependent  pulse  train.  The  algorithm  used  in  the 
software  combines  the  spatial  configuration  of  the  beam  with  the  temporal  configuration. 

The  other  important  parameter  in  the  pump  algorithm  is  the  input  power.  Care  must  be  taken  in 
determining  the  power  per  pulse  in  a  time-dependent  algorithm.  The  software  itself  ensures  that 
the  correct  amount  of  power  is  available  to  each  time  increment  in  the  run.  That  is,  if  the  power 
input  is  set  to  one  watt  and  the  minimum  time  increment  is  1  ps,  then  the  software  provides  1 
pW  to  each  time  increment.  Once  the  pump  algorithm  acts,  this  power  will  either  be  put  into  the 
material  or  not,  thus  the  peak  power  is  always  the  same  but  the  average  power  will  depend  on 
duty  cycle.  This  is  the  analog  of  a  chopped  pump  beam  where  a  constant  amount  of  CW  pump 
power  is  chopped  to  achieve  a  quasi-CW  modulated  pump  beam. 

If  the  goal  is  to  model  a  Q-switched  pump  the  average  power  is  generally  the  same  (as  long  as 
the  reciprocal  of  the  pulse  repetition  frequency,  1/PRF,  is  less  than  the  relaxation  lifetime  of  the 
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pumping  laser)  but  the  peak  pump  power  depends  on  PRF.  Of  note  in  this  case  is  the  fact  that 
most  Q-switched  pump  lasers  used  to  pump  Cr2+  lasers  have  a  pulse  width  of  100  ns  or  less 
which  is  much  smaller  than  the  smallest  practical  time  increment  in  our  modeling  program.  We 
have  approximated  Q-switched  pumping  in  our  model  by  keeping  average  power  constant,  using 
a  constant  low  duty  cycle,  and  varying  PRF. 

3. 1.2. 3  Minimum  time  increment 

The  FemLab®  time  dependent  solver  attempts  to  find  a  solution  at  discrete  time  increments 
throughout  the  total  time  period.  The  time  increment  chosen  by  the  user  should  strike  a  balance 
between  model  accuracy  and  solution  time  required.  It  should  also  fully  resolve  the  pulses  in  the 
pump  algorithm.  Time  increment  variations  of  10,  5,  2  and  1  ps  were  modeled  for  a  10-kHz  PRF 
and  15%  duty  cycle  in  order  to  find  the  optimal  time  increment  setting.  Figure  5,  below,  shows 
the  results  of  the  z-average  temperature  at  the  x-y  center  point  for  the  different  time  increments. 


Figure  5:  The  z-average  temperature  at  the  x-y  center  point  for  the  different  time  increments  over  a  total  of  4 
ms 


The  difference  between  1  ps  and  2  ps  is  about  5.5%  of  the  average  temperature  rise  (above  293 
K).  To  study  the  differences  between  1  and  2  ps,  these  models  were  run  for  8  ms. 
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Figure  6:  The  z-average  temperature  at  the  x-y  center  point  comparison  for  two  different  time  increments  for 

a  total  of  8  ms 


As  can  be  seen  in  Figure  6,  the  1  gs  model  seems  to  be  the  better  model  due  to  fluctuations  in  the 
2  gs  model,  although  both  are  converging  well  at  the  end  of  8  ms,  to  within  2%  of  the  average 
temperature  rise.  The  smallest  time  increment  feasible  due  to  time  and  memory  constraints  was 
1  gs.  Modeling  results  for  larger  time  increments  were  inconsistent  throughout  the  model  runs 
but  converged  on  the  1  gs  results.  This  indicated  that  1  gs  was  the  optimal  time  increment  to  use 
for  all  models.  This  would  indicate  that  the  conductivity  of  ZnSe  causes  any  changes  in  the  heat 
equation  to  take  about  1  gs  to  propagate  through  the  crystal  enough  to  effect  the  finite  element 
analysis.  This  could  possibly  change  then  by  changing  the  number  of  elements  (thus  the  size  of 
each  element). 

Due  to  time  and  memory  constraints,  it  is  difficult  to  run  a  1-gs  time  increment  model  for  8  ms 
or  longer  total  time.  It  would  be  useful  to  be  able  to  model  longer  total  time  periods  by  modeling 
shorter,  discrete  time  periods.  The  last  solution  of  the  previous  run  would  be  used  as  the  initial 
condition  for  the  next  model  with  the  results  merged  into  a  continuous  solution.  The 
comparison,  in  Figure  7,  of  a  continuous  8-ms  model  to  two  discrete  4-ms  models  showed  that 
this  method  is  acceptable  as  the  differences  in  results  between  the  two  methods  were  negligible 
in  comparison  to  the  modeled  values. 
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Figure  7:  Comparison  of  single  continuous  model  run  to  combine  discrete  runs  showing  acceptability  of  using 
discrete  model  runs  for  computational  convenience 


The  1-jj.s  time  increment  model  was  also  run  for  16  ms  to  see  how  closely  the  temperatures 
approached  an  asymptote.  This  data  can  be  seen  in  Figure  8. 


Figure  8:  16  ms  modeling  run  at  1  ps  increments  to  determine  asymptotic  behavior 

At  8  ms,  the  z-average  temperature  rise  (defined  as  the  average  temperature  rise  over  a  line 
through  the  depth  of  the  sample  in  the  propagation  (z)  direction)  is  within  0.2  K  of  the  value  at 
16  ms.  This  difference  is  less  than  5%  of  the  z-average  temperature  rise.  This  number  is  also 
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within  5%  of  the  steady  state  average  temperature  rise  of  9.3  K.  This  confirmed  that  the 
maximum  total  time  needed  to  run  the  models  to  achieve  asymptotic  behavior  is  8  ms. 

3. 1.2. 4  Determination  of  ancillary  parameter  values 

The  goal  of  the  time  dependent  modeling  is  to  study  the  effects  of  different  pulse  repetition 
frequencies  (PRFs)  and  duty  cycle  on  both  the  magnitude  of  thennal  lensing  and  the  overall 
average  temperature  rise  in  the  laser  material.  There  are  a  large  number  of  parameters  in  models 
such  as  this  (geometry,  time  increment,  mesh  size,  etc)  and  all  of  them  except  for  PRF  and  duty 
cycle  must  be  kept  constant  from  run  to  run  in  order  to  make  valid  comparisons.  A  large  portion 
of  the  modeling  work  involved  determination  of  appropriate  values  for  these  ancillary 
parameters.  The  optimal  value  is  one  which  produces  accurate  and  precise  results  while  still 
allowing  for  reasonable  solution  times. 

In  previous  work,  the  maximum  portion  of  the  laser  material  geometry  that  needed  to  be  included 
in  the  simulation  was  detennined  [2].  It  was  postulated  that  areas  of  the  model  farther  from  the 
center  than  approximately  twice  the  pump  beam  radius  would  be  irrelevant  to  the  model  solution 
and  could  thus  be  excluded.  A  set  of  basic  face  cooled  models  (modeled  as  parallelepipeds  for 
simplicity)  was  constructed  where  the  x  and  y  dimensions  varied  from  2  mm  to  8  mm  while  the 
pump  beam  w  (see  eq.  (3.1))  was  0.4  mm.  The  thennal  lensing  power  was  found  to  not  vary 
significantly  in  this  range  with  size  and  4-mm  sides  produced  an  upper  bound  on  thennal  lensing 
power  and  maximum  average  temperature  rise.  From  these  models,  it  was  determined  that  the  4- 
mm  x  andy  size  was  optimal. 

Due  to  limited  computer  memory  and  CPU  resources,  it  was  desirable  to  run  the  models  with  as 
few  degrees  of  freedom  (a  measure  of  the  model  complexity)  as  possible  without  sacrificing 
model  quality.  To  insure  that  model  complexity  was  not  lowered  to  the  point  of  affecting  the 
results,  the  standard  4-mm  model  was  run  at  varying  degrees  of  freedom.  The  results  showed 
that  the  thennal  lensing  power  and  maximum  average  temperature  rise  leveled  off  quickly  and 
that  as  long  as  some  reasonable  mesh  size  was  chosen,  the  results  would  be  acceptable.  With 
this  software  package,  solution  time  was  not  the  main  impediment,  memory  consumption  was. 
The  most  complex  model  that  available  1  GB  RAM  computers  could  run  was  about  220000 
degrees  of  freedom  due  to  memory  allocation  issues. 

Once  the  appropriate  x-y  dimensions  and  model  mesh  complexity  were  established,  the  geometry 
of  the  sample  itself  was  studied.  Due  to  symmetry  properties,  it  was  beneficial  to  use  thin 
cylinders  instead  of  parallelepipeds  in  the  models.  A  4-mm  diameter  cylinder  was  compared  to 
the  4-mm  by  4-mm  parallelepiped  to  gauge  any  differences.  Modeling  showed  that  differences 
in  thennal  lensing  power  and  maximum  average  temperature  rise  between  the  disk  and 
parallelepiped  were  two  orders  of  magnitude  smaller  than  the  lensing  power  itself.  Thus  the 
difference  between  the  disk  and  parallelepiped  geometries  was  negligible  when  looking  at 
thermal  lensing  power  and  maximum  average  temperature  rise. 

3.2  Model  Results  and  Analysis 

Once  the  geometric  models  were  constructed  and  configured  as  required,  the  finite  element 
solver  was  initiated.  FemLab®  provided  the  solutions  to  the  each  time  increment  as  data 
structures  which  contained  all  of  the  relevant  temperature  and  geometry  data.  It  also  provided  a 
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visual  display  of  the  solution  in  its  post-processing  mode.  These  solutions  showed  the 
temperature  distribution  along  a  slice  along  one  axis  through  the  origin.  As  an  example,  the 
solution  for  the  4-mm  basic  face  cooled  thin  disk  design  can  be  seen  in  Figure  9  and  Figure  10. 


Figure  9:  Simulation  solution  of  4  mm  diameter,  1  mm  thick  face  cooled  thin  disk  pumped  with  a  Gaussian  (n 
=  1)  beam  with  a  =  0.4  mm.  The  temperature  ranges  from  293  (dark)  to  306  (light). 


Figure  10:  Temperature  gradient  surface  for  previous  4  mm  face  cooled  thin  disk  solution 

After  a  temperature  solution  for  each  mesh  point  was  determined  for  every  time  increment,  the 
raw  temperature  data  needed  to  be  processed  into  a  more  useful  form.  In  this  research,  we  were 
interested  in  two  types  of  temperature  information.  The  maximum  z-average  temperature  rise  in 
the  laser  material  has  an  effect  on  the  fluorescence  lifetime  while  the  radial  gradient  of  the  z- 
average  temperature  is  an  indication  of  the  magnitude  of  the  thermal  lensing.  To  obtain  the  z- 
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average  temperature,  the  data  for  each  time  increment  must  be  exported  and  then  z-averaged. 
The  radial  temperature  gradient  is  important  because  of  the  high  thenno-optic  coefficient  of  the 
laser  material.  As  temperature  increases,  so  does  the  refractive  index  which  means  that  a 
temperature  gradient  will  produce  an  effective  quadratic  thennal  lens.  The  power  of  the  thermal 
lens  is  related  to  how  sharply  the  gradient  peaks.  The  pump  cycle  has  several  interesting  times  at 
which  the  thermal  gradient  could  be  measured.  These  would  be  mid-way  through  the  pulse,  at 
the  peak  of  the  pulse,  mid-way  through  the  rest  portion  of  the  cycle  and  just  prior  to  the  pulse 
starting.  Because  lasing  will  start  at  the  beginning  of  the  pump  pulse,  that  is  where  we  measured 
the  radial  temperature  gradient. 

As  mentioned  in  a  previous  section,  most  pulsed  lasers  fall  into  one  of  two  categories.  Lasers  we 
will  refer  to  as  chopped  lasers  have  a  fixed  peak  power  and  an  average  power  which  varies  with 
duty  cycle,  all  of  which  are  mostly  independent  of  PRF.  This  would  be  analogous  to  directly 
modulating  or  physically  chopping  a  CW  fiber  laser  such  as  an  IPG  TmiFiber  laser  used  in  this 
research  as  a  pump  source.  Lasers  we  will  refer  to  as  Q-switched  lasers  have  a  fixed  average 
power,  a  fixed  pulse  width  and  a  peak  power  and  duty  cycle  which  vary  with  PRF.  This  would 
be  represented  by  the  Tm:Ho:YLF  laser  used  in  previous  works  as  a  pump  source  for  gain 
switched  lasing  of  Cr  :ZnSe.  A  final  category  of  laser  which  can  be  modeled  is  a  hybrid  system 
where  the  average  power  is  fixed  as  in  the  Q-switched  case,  but  the  pulse  width  varies  with  the 
duty  cycle  as  in  the  chopped  case.  Since  no  experimental  analog  to  this  was  in  use,  it  was  not 
modeled.  In  our  modeling  calculations  below  we  assumed  1  W  of  cw  pump  power;  other  powers 
can  be  considered  by  remembering  that  temperature  rise  and  thennal  lensing  power  scale  linearly 
with  pump  power.  We  also  assumed  a  1  -mm  thick  disk  4  mm  in  diameter  which  absorbed  all 
pump  power. 

3.2.1  Chopped  lasers 

In  a  chopped  laser,  the  peak  power  is  fixed,  the  average  power  is  determined  by  the  duty  cycle 
and  both  are  independent  of  the  PRF.  CW  operation  corresponds  to  1  W  peak  and  average 
power  while  25%  duty  cycle,  for  example,  is  1  W  peak  power  and  0.25  W  average  power.  For  a 
1-kHz  PRF  at  varying  duty  cycles,  the  z-average  temperature  at  the  center  point  can  be  seen  in 
Figure  1 1 .  Note  that  the  temperature  rise  does  not  completely  return  to  the  baseline  between 
pulses  but  the  average  temperature  does  reach  an  asymptotic  value  after  about  8  ms  of  pumping. 
All  duty  cycles  other  than  cw  operation  show  a  saw-tooth  (l-exp(-t/x))  rise  and  fcxpf-t/x))  decay 
at  all  times  where  x  is  the  relaxation  lifetime.  Also  note  that  just  1  W  cw  can  result  in  a  9  K 
temperature  increase;  multiple-watt  pumping  could  easily  reach  temperatures  where  nonradiative 
relaxation  would  become  severe. 
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Figure  11:  Center  point  z-average  temperature  for  1  kHz  PRF  at  varying  duty  cycles. 


The  radial  thermal  profiles  for  1  kHz  PRF  and  varying  duty  cycles  are  shown  in  Figure  12. 
Table  3  shows  calculations  of  thermal  lensing  power  for  the  different  duty  cycles.  Obviously, 
thermal  lensing  should  increase  with  duty  cycle  (more  average  power);  in  fact,  it  increases 
quadratically  with  duty  cycle  due  to  increased  pump  power  and  decreased  time  to  cool. 
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Figure  12:  Radial  temperature  gradient  for  1  kHz  PRF  at  varying  duty  cycles 


Table  3:  Thermal  Lensing  Power  at  1  kHz  PRF,  1  W  peak  power,  Chopped  Laser 


Duty  cycle  (%) 

Thermal  lens  power  (m'1) 

10 

0.23 

25 

0.63 

35 

0.88 

50 

1.55 

75 

3.06 

100 

6.00 

These  modeling  results  predict  that  lower  duty  cycles  should  better  support  lasing  through  lower 
temperatures  and  lower  thermal  lensing.  This  is  supported  experimentally  as  seen  by  the  results 
in  Figure  13  [37].  At  low  duty  cycles  the  laser  output  power  slope  efficiency  remained  linear  for 
all  average  pump  powers  up  to  5  watts.  However,  the  high  duty  cycle  results  showed  rollover, 
i.e.  reduced  slope  efficiency  above  4.5  W  average  pump  power  for  50%  duty  cycle  and  above  3.7 
W  for  100%  duty  cycle  operation.  The  threshold  for  rollover  is  not  simply  a  function  of  average 
power  because  the  pump-off  times  allow  the  laser  material  temperatures  to  cool  below  the  cw 
average.  The  decrease  in  slope  efficiency  with  increasing  duty  cycle  is  likely  due  to  the 
increased  average  temperature  of  the  sample  thus  reducing  Cr2+  energy  storage  time.  The  onset 
of  rollover  is  likely  to  be  due  to  onset  of  resonator  instability  from  thermal  lensing. 
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Figure  13:  Experimental  results  for  1  kHz  PRF  chopped  laser  with  varying  duty  cycles.  Pump  power  is 
average  power. 


Figure  14:  Experimental  results  for  fixed  PRF  and  duty  cycles  for  varying  heat  sink  temperatures 

The  model  predicts  that  changing  the  heat  sink  temperature  should  only  effect  the  slope 
efficiency  and  not  the  power  drop-off  associated  with  thermal  lensing  or  the  threshold  power.  In 
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order  to  verify  the  influence  of  baseline  heat  sink  temperature  on  the  laser  performance,  a  series 
of  slope  efficiency  measurements  were  made  at  different  temperatures.  The  other  parameters 
were  chosen  based  on  previous  measurements  that  indicated  thermal  lensing  at  higher  powers. 
As  Figure  14  indicates,  the  power  drop-off  at  higher  pump  powers  due  to  thermal  lensing 
remains  similar  between  the  run  as  does  the  threshold  power.  The  slope  efficiency  decreases  as 
the  heat  sink  temperature  increases.  Heat  sink  temperatures  below  15°  C  caused  condensation 
and  were  not  suitable  for  laser  use  without  dry  air  purging. 

3.2.2  Q-switched  Lasers 

In  a  system  where  the  average  power  and  duty  cycle  are  fixed,  variations  in  PRF  have  little  effect 
on  the  overall  time-average  temperature  at  the  center  of  the  material.  This  is  because,  over  a 
period  of  time  which  is  large  with  respect  to  the  pulse  period,  the  same  amount  of  heat  is  being 
put  into  the  laser  with  the  same  relative  amount  of  rest  time  (due  to  the  fixed  duty  cycle).  Figure 
15  shows  the  z-average  temperature  for  15%  duty  cycle  pumping  and  varying  repetition  rate. 
The  model  implicitly  assumes  instantaneous  heating  with  absorption  of  pump  power.  It  shows 
that  during  the  pumping  time  (which  is  when  gain  induced  Cr“  lasing  would  be  occurring),  the 
average  temperature  goes  higher  as  the  repetition  rate  decreases.  This  contradicts  the 
experimental  observation  that  average  Cr"  output  power  is  constant  with  increasing  PRF  until  it 
approaches  10  kHz  and  then  output  declines.  However,  pumped  Cr2+  ions  do  not  immediately 
transfer  their  stored  energy  to  the  crystal  lattice;  either  stimulated  emission  or  transition  back  to 
the  ground  state  must  occur.  Thus  the  temperatures  at  the  beginning  of  the  pumping  are  likely  to 
be  more  representative  of  the  temperature  during  lasing  for  the  Q-switched  case.  In  that  case,  the 
minimum  temperature  increases  with  PRF  reaching  the  cw  limit  as  PRF  goes  to  infinity. 

The  radial  temperature  gradients  were  extracted  at  the  low  temperature  points  of  each  PRF  cycle 
in  Figure  15  at  times  when  the  average  temperature  was  close  to  steady  state.  This  is  based  on 
the  assumption  that  thermal  lensing  effects  were  most  important  at  the  beginning  of  the  pump 
cycle.  These  gradients  can  be  seen  in  Figure  16  and  the  calculated  thermal  lens  powers 
associated  with  them  are  listed  in  Table  4. 


18 


Temperature  (K)  2  Temperature  (K) 


Time(s)  x1Q-3 

15:  Center  point  z-average  temperature  for  15%  duty  cycle  with  varying  PRF 


Figure  16:  Radial  temperature  gradient  for  15%  duty  cycle  at  varying  PRFs. 
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Table  4:  Thermal  Lensing  Power  at  1  kHz  PRF,  1  W  peak  power,  Q-switched  Laser 


PRF  (kHz) 

Thennal  lens  power  (m"‘) 

0.5 

1.40 

1 

2.43 

2 

3.60 

5 

4.84 

10 

5.37 

15 

5.59 

20 

5.97 

3.3  Modeling  Conclusions 

Time-dependent  modeling  of  thermal  effects  in  laser  materials  is  a  viable  option  for 
understanding  the  influences  of  laser  parameters  such  as  PRF  and  duty  cycle  on  laser 
perfonnance.  The  model  is  in  agreement  with  experimental  results  for  the  chopped  pumping 
case  but  does  not  directly  agree  with  Q-switched  pumping  because  the  time  delay  between 
absorption  and  heat  transfer  to  the  host  material  is  not  accounted  for  in  our  model.  Future 
modeling  work  will  focus  on  clarifications  of  the  chopped  laser  results  and  on  development  of  a 
way  to  model  the  time  delay  between  Q-switched  pumping  and  heating  of  the  host  material. 
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Chapter  4 

Chromium-doped  Zinc  Selenide  Fiber  Laser 

4.1  Introduction 

As  previously  mentioned,  a  number  of  military  and  scientific  applications  require  laser  sources  in 
the  mid-IR  region.  Due  to  the  high  stress  environment  in  which  these  types  of  laser  sources  are 
likely  to  be  used,  they  must  be  compact,  rugged  and  solid-state.  Vibration,  dust,  limited  power 
and  little  to  no  cooling  are  all  issues  that  must  be  addressed  in  any  architecture  of  this  type. 
Doped  fiber  amplifier  laser  systems  are  optimal  candidates  for  this  type  of  use  due  to  their 
inherent  properties.  In  addition  to  environmental  ruggedization,  multiple  fibers  could  be  used 
with  a  single  pump  source  or  vice-versa.  This  ability  to  possibly  distribute  laser  power  to 
various  points  around  an  aircraft  body  using  a  central  pump  source  could  be  extremely  useful  in 
applications  which  require  multiple  output  ports  such  as  IRCM  systems. 

Silica  fiber  has  long  been  used  in  the  telecom  industry  as  both  a  transport  medium  and  for  all- 
optical  amplification  applications  such  as  Erbium  Doped  Fiber  Amplifiers  (EDFAs).  An  EDFA 
is  essentially  a  normal  optical  fiber  waveguide  that  has  been  lightly  doped  with  a  rare-earth 
element,  namely  erbium  (possibly  co-doped  with  ytterbium  to  increase  the  pumping  efficiency 
and  gain).  In  a  typical  EDFA,  980  mn  photons  are  used  to  create  1550  nm  photons  which  serve 
to  amplify  signal  photons  at  the  same  wavelength. 

Using  this  same  idea,  silica  fibers  are  useful  for  the  construction  of  fiber  lasers.  These  glasses 
have  a  high  damage  threshold,  their  surface-to-volume  ratio  is  high  leading  to  easy  heat 
dissipation,  the  gain  medium  is  nicely  contained  in  a  waveguide  and  current  fiber  manufacturing 
capabilities  make  available  a  wide  range  of  quality  fibers.  Due  to  the  nature  of  silica  (with 
strong  absorption  starting  around  2  pm),  silica  fiber  lasers  have  generally  been  limited  to  sub-2 
pm  applications.  Commercially  available  lasers  are  commonly  found  using  Er  (1.5  pm)  and  Tm 
(1.9  pm). 

Expanding  on  the  available  methods  for  thennal  lensing  mitigation,  one  goal  was  to  demonstrate 
a  chromium-doped  zinc  selenide  (Cr~  :ZnSe)  fiber  amplifier.  One  method  demonstrated  in  the 
literature  for  fabricating  optical  fibers  from  crystalline  or  polycrystalline  laser  hosts  was  to  using 
a  standard  fiber  drawing-tower  method  with  a  larger  diameter,  pre-fiber  form  or  preform  which 
was  constructed  using  the  rod-in-tube  (RIT)  method  [49,  50].  Using  this  method,  a  silica  rod 
was  tapered  at  one  end  and  had  a  hole  drilled  in  the  other.  A  Cr2  :ZnSe  crystal  rod  doped  to  the 
appropriate  level  was  inserted  in  this  hole,  which  was  then  capped  with  silica  to  contain  the 
crystal.  The  dimensions  of  the  preform  and  crystal  were  designed  to  give  the  desired 
core/cladding  size  ratio  after  drawing.  The  melting  point  of  zinc  selenide  is  1526°C  while  the 
softening  point  of  silica  is  1670°C  which  should  have  allowed  for  consistent  drawing.  The 
coefficient  of  thennal  expansion  of  zinc  selenide  is  much  higher  than  that  of  silica,  so  heating  of 
the  prefonn  was  done  slowly  to  avoid  thennal  shock.  Further,  it  was  hoped  that  the  preform 
could  be  manufactured  to  produce  modified  fiber  designs  such  as  a  dual-core  design  to  allow 
pumping  in  the  outer  core  with  signal  amplification  and  guiding  in  the  inner  core. 
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4.2  Experimental  Results 

Discussions  were  initiated  with  researchers  at  Clemson  University  to  fabricate  a  prefonn  like  the 
one  described  above.  In  addition,  Cr2+:ZnSe  with  the  proper  dimensions  to  fit  in  the  perfonn 
hole  was  also  needed.  Previous  samples  of  chromium  zinc  selenide  were  purchased  from  a 
commercial  vendor  which  had  subsequently  gone  out  of  business,  so  a  search  for  a  new  vendor 
was  initiated.  One  of  our  other  sources  of  Cr  :ZnSe  in  the  past  was  MIT/Lincoln  Labs.  Costs 
for  resumption  of  chromium  doping  by  this  group  was  discussed  and  fell  in  the  $65K  range 
which  was  considered  too  high  for  the  risk  level  of  this  project.  After  discussions  with  Sergey 
Mirov’s  group  at  the  University  of  Alabama  Birmingham  (UAB),  it  was  decided  that  they  could 
provide  the  types  of  samples  needed. 

Prior  to  the  procurement  of  new  Cr“  :ZnSe,  two  uncoated,  chromium-doped  zinc  selenide 
samples  of  dimensions  10x9x2  mm  and  8x10x2  mm  respectively  were  sent  to  Clemson  for  use  in 
fiber  pulling  experiments.  These  samples  were  crushed  up  and  used  as  core-filler  for  a  fiber 
prefonn  (see  Figure  17).  These  were  to  be  used  in  place  of  the  rods  in  order  to  begin  producing 
fiber  for  testing.  These  samples  had  volumes  of  180  mm  and  160  mm  respectively.  The 
planned  50  mm  long  x  5  mm  diameter  core  rod  would  have  a  volume  of  about  982  mm3, 
considerably  larger  than  the  given  samples.  In  order  to  get  a  fiber  draw  of  an  appreciable  length 
and  good  quality,  both  samples  were  used  for  one  draw.  The  goal  was  to  get  a  125  pm  diameter, 
25  pm  core  fiber.  The  fiber  prefonn  used  was  an  F300  tube,  6  mm  ID  x  30  mm  OD  x  400  mm 
length.  Over  500  m  of  125  pm  diameter  fiber  was  produced  from  the  preform.  This  draw  was 
done  under  vacuum  to  minimize  air  inclusion. 


^  Vacuum  sealed 

Figure  17:  Fiber  preform  design 

Initial  analysis  indicated  that  the  draw  was  not  successful  as  most  of  the  Cr  :ZnSc  seemed  to 
have  volatilized  and  moved  further  up  the  preform  as  heat  was  applied  during  drawing.  As  this 
material  reached  a  cooler  prefonn  zone,  it  crystallized  as  seen  in  Figure  18.  No  Cr2+  seems  to  be 
present  in  the  draw-end. 

Because  of  this,  subsequent  draws  were  re-designed  to  fill  any  open  portion  of  the  preform  not 
occupied  by  laser  host  with  a  silica  rod  and  seal  it. 


V  2+ 

^  Crushed  Cr  :ZnSe 
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rystallized  Cr2+:ZnSe 


Figure  18:  Preform  remaining  after  first  fiber  draw  attempt 

Four  of  the  six  cylindrical  (~5mm  diameter  by  48  mm  length)  Cr  :ZnSe  samples  from  UAB 
were  received  and  were  taken  to  Clemson  University’s  Center  for  Optical  Materials  Science  and 
Engineering  Technologies  (COMSET)  for  fiber  pulling.  The  second  fiber  draw  attempt  began 
with  preform  construction  using  COMSET’s  cutting  edge  lathe  system.  As  seen  in  Figure  19,  a 
Heraeus  quartz  tube,  30  mm  OD,  6  mm  ID  and  about  300  mm  long  was  first  heated  on  one  end 
while  the  ID  was  enlarged  slightly  using  a  carbon  reamer.  This  was  done  to  achieve  a  ‘funnel’  to 
facilitate  the  insertion  of  materials  into  the  preform.  This  was  then  attached  to  a  long,  thin- 
walled  hollow  handle  tube  for  the  chuck  (‘C’)  end.  The  purpose  of  this  piece  was  to  provide 
distance  between  the  furnace  and  the  chuck  assembly  in  the  drawing  tower  as  well  as  to  allow  for 
connection  to  a  vacuum  manifold.  A  solid  piece  of  quartz,  6  mm  OD  by  —20  mm  long  was  then 
inserted  into  the  furnace  (‘F’)  end  of  the  Heraeus  tube  and  sealed  by  attaching  a  solid  20  mm  OD 
by  -200  mm  long  rod. 


Figure  19:  Fiber  preform  assembly,  step  one 

The  solid  furnace-end  rod  was  then  melted  near  the  tube  and  pulled  to  achieve  a  tapered  end. 
The  Cr  :ZnSe  cylinder  was  then  slid  down  the  hollow  chuck  end,  through  the  funnel  and  down 
into  place  at  the  sealed  furnace  end.  A  6  mm  OD  by  -400  mm  long  ‘tamper’  tube  was  then 
inserted  behind  the  Cr2  :ZnSe  sample  to  reduce  volatilization.  The  preform  was  evacuated  and 
sealed  at  the  very  chuck-end  of  the  hollow  tube.  The  final  preform  can  be  seen  in  Figure  20. 
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Figure  20:  Preform  construction,  step  two 

The  preform  was  lowered  into  the  drawing  tower  furnace  zone  and  a  typical  drawing  routine  was 
followed.  The  prefonn  subsequently  exploded  while  outside  the  draw  zone  which  was  at  1700° 
C.  The  prefonn  had  Tight  bulbed’  or  blown  out  into  the  furnace  zone.  The  resulting  damage  to 
the  prefonn  can  be  seen  in  Figure  21.  Due  to  the  nature  of  zinc  selenide,  an  emergency  stop  was 
perfonned  and  the  room  evacuated  until  the  furnace  completely  cooled.  An  entire  day  of  cleanup 
as  well  as  furnace  element  and  heat  shield  replacement  followed.  Zinc  selenide  begins 
sublimating  at  temperatures  below  the  “drop  point”  (point  at  which  the  fiber  can  start  being 
pulled)  [51]  and  that  this  is  the  source  of  the  pressure  which  caused  the  explosion.  Further  fiber 
draws  were  placed  on  hold  until  a  safer  draw  method  could  be  devised. 


Figure  21:  Exploded  preform  from  second  fiber  draw  attempt 
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Fluorescence  decay  lifetime  measurements  were  made  on  the  first  fiber  prefonn  to  look  for  signs 
of  Cr2+  ions.  In  order  to  access  the  desired  areas  for  analysis,  the  preform  was  roughly  polished 
to  remove  oxidized  silica  and  cleaved  near  the  crystalline  material,  as  seen  in  Figure  22.  The 
~15  mJ  excitation  pulses,  generated  by  a  1  Hz,  flashlamp-pumped,  AO  Q-switched 
Tm:Ho:Cr:YAG  2.1  pm  laser,  were  aimed  at  both  the  crystalline  material  location  as  well  as  a 
clear  quartz  location  for  a  control.  The  lifetime  results  can  be  seen  in  Figure  23.  The  pulse 
shapes,  as  seen  by  an  extended  range  InGaAs  PIN  diode  showed  that  the  core  had  a  noticeable 
decay  after  the  pulse  when  compared  to  the  control.  Using  a  MatLab  curve  fitting  algorithm,  the 
lifetime  of  the  core  was  found  to  be  3.8  ps,  which  is  very  much  in  line  with  the  measured 
lifetime  of  Cr2+  in  other  samples  from  the  same  doping  run.  This  indicated  a  high  likelihood  that 
Cr“  :ZnSe  can  melt  and  re-crystallize  while  keeping  the  chromium  in  the  important  2+  valence 
state  required  for  lasing  in  the  2-3  pm  region. 


Figure  22:  Polished  and  cleaved  preform,  ready  for  spectroscopy 
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Figure  23a:  Decay  pulses  from  preform;  7b:  Lifetime  measurement  of  core 

Spectroscopy  was  performed  on  the  first  fiber  preform  using  an  Acton  SpectraPro  150 
monochromator,  SpectraHub  detector  interface,  an  extended  range  InGaAs  detector  and 
SpectraSense  software.  The  results,  seen  in  Figure  24,  show  definitively  that  the  material  which 
melted  and  re-solidified  is  indeed  Cr  :ZnSe.  The  response  of  a  mix  of  core  and  cladding  show 
the  expected  strong  absorption  of  the  silica  in  the  range  above  2200  nm. 
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Figure  24:  Spectroscopy  results  for  different  postions  on  fiber  preform 

A  spectroscopy  system,  seen  in  Figure  25,  for  use  with  the  experimental  chromium  fibers  was 
assembled.  The  system  was  tested  for  calibration  at  1900  nm  using  an  IPG  thulium  liber  laser 
and  three  older  samples  of  Cr“  :ZnSe.  The  pump  beam  was  chopped  with  a  SciTech  optical 
chopper  running  at  76  Hz  and  the  detector,  an  Infrared  Associates  InSb,  was  read  using  a 
Stanford  Research  Systems  SR830  lock-in  amplifier  connected  to  the  SpectraHub.  The  slit 
widths  were  set  at  1.17  pm  while  the  pump  power  was  less  than  30  mW.  Data  taken  with  this 
setup  confirmed  that  no  detectable  Cr  :ZnSe  was  contained  in  the  liber  drawn  from  the  first 
preform.  This  confirms  the  speculation  above  that  any  Cr2+:ZnSe  was  vaporized  and  moved  up 
and  out  of  the  draw  zone  before  it  could  be  incorporated  into  the  fiber  core. 
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Figure  25:  Fiber  spectroscopy  setup 

COMSET  proposed  a  new  fiber  prefonn  construction  and  drawing  experiment  based  on  technical 
discussions  of  the  volatilization  problem.  Two  new  preforms  were  constructed,  one  using  a 
much  shorter  solid  Cr2  :ZnSc  cylinder  (20  mm  or  less)  and  the  other  using  a  smaller  diameter  of 
crushed  Cr  :ZnSe,  as  seen  in  Figure  26.  Using  a  shorter  cylinder  should  allow  the  preform  to 
have  all  of  the  ZnSe  in  the  hot  zone  of  the  furnace  at  the  same  time  thus  preventing  sublimation 
and  re-crystallization  at  another  location  along  the  ZnSe  rod.  The  re-crystallization  was  to  seal 
the  vapor  into  a  small  volume.  Reducing  the  volume  of  ZnSe  being  heated  should  also  reduce  the 
vapor  pressure.  The  crushed,  small  diameter  sample  was  to  have  low  vapor  pressure  due  to  low 
volume  and  the  crushed  nature  of  the  material  should  allow  vapor  to  escape  the  hot  zone. 
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Figure  26:  New  preform  designs  a)  solid,  short  core  and  b)  long,  small  diameter  powdered 
core 

New  fiber  draws  incorporating  the  new  designs  were  attempted.  As  with  all  previous  attempts 
save  the  first  one,  the  prefonn  exploded  prior  to  the  drop  point. 
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4.3  Conclusions 

Further  fiber  draws  of  Cr“  :ZnSe  cores  were  not  attempted  due  to  the  cost  and  hazards  associated 
with  preform  explosions.  It  was  hypothesized  that  the  first  draw  only  worked  because  the  small 
amount  of  chromium  zinc  selenide  powder  in  the  first  prefonn  was  completely  volatilized  up 
onto  the  walls  of  the  upper  cold  section  prior  to  drop,  leaving  no  material  to  explode  and 
resulting  in  no  Cr“  :ZnSe  incorporation  into  the  fiber.  Future  work  should  shift  to  other 
techniques.  Penn  State  University  has  shown  promising  developments  in  depositing  metals  and 
semiconductors  in  hollow  core  fibers  [52].  It  may  also  be  possible  to  fabricate  waveguide 
structures  in  bulk  Cr2+:ZnSe  using  ultrashort  laser  pulses  [53]. 
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Chapter  5 

High-power,  Widely-tunable  Cr2+:ZnSe  Lasers 

5.1  Introduction 

This  chapter  presents  work  undertaken  to  implement  lessons  learned  through  modeling  and 
previous  work  [2,  3,  33,  37,  47,  54]  with  face-cooled  Cr2+:ZnSe  lasers  to  construct  high-power, 
widely-tunable  laser  systems.  The  previous  work  illustrated  that  the  thermal  issues  limiting 
Cr"  :ZnSe  laser  performance  can  be  addressed  through  appropriate  choices  of  crystal  geometry 
(helping  to  reduce  thermal  gradients),  a  longer  wavelength  pump  source  (reducing  quantum- 
defect)  and  improved  heat  removal;  all  of  which  help  to  reduce  the  thermal  load  in  the  gain 
element.  However,  even  when  these  steps  are  taken  to  mitigate  the  thermal  problems,  the  lasers 
are  often  still  limited  in  their  output  characteristics.  One  solution  to  this  problem  is  to  use  a 
master  oscillator  power  amplifier  (MOP A)  architecture  where  the  MO  design  is  chosen  to 
emphasize  a  desirable  operating  characteristic  such  as  beam  quality,  linewidth  or  tunability  while 
the  amplifier  section  accomplishes  the  power  scaling  without  degrading  performance. 

This  research  investigated  MOPA  laser  perfonnance  using  polycrystalline  Cr"  :ZnSe  gain  media 
in  the  conventional  slab  geometry  for  both  the  oscillator  and  the  amplifier.  We  developed  an 
improved  amplifier  gain  model  which  incorporates  variations  in  pump  and  signal  beam 
parameters.  We  investigated  a  variety  of  MO  configurations  and  their  performance  as  part  of  a 
versatile  MOPA-configuration  laser  system.  We  demonstrated  greater  than  10  W  of  pure 
continuous  wave  (CW)  output  and  compared  the  system  performance  to  the  theoretical  model. 
System  versatility  was  demonstrated  by  multi-watt,  narrow  linewidth  output  tunable  over  a  range 
greater  than  400  mn.  Beam  quality  and  spectral  content  were  also  investigated. 

5.2  Amplifier  Model 

In  order  to  design  a  MOPA  Cr"  :ZnSe  laser,  a  theoretical  model  was  developed  and  used  to 
detennine  possible  amplifier  performance  for  a  given  set  of  physical  parameters.  We  considered 
optical  pumping  co-propagating  and  co-aligned  with  an  input  signal  beam.  Absorption  by  the 
Cr“  ions  is  proportional  to  pump  intensity  so  both  pump  intensity  and  absorbed  power  will 
decrease  along  the  propagation  direction.  In  addition,  both  the  pump  and  signal  are  expected  to 
have  some  non-uniform  transverse  distribution  which  we  will  model  as  well. 

Any  amplifier  model  must  be  based  on  known  or  measured  material  parameters  involved  in 
absorption  and  emission.  The  absorption  and  emission  cross-sections  for  Cr2  :ZnSe  vary  in  the 
literature  (for  example,  between  90  x  10'  and  130  x  10'  cm"  for  peak  emission  cross-section). 
This  model  will  use  values  published  by  DeLoach  et.  al.  [27]  but  also  acceptable  would  be  those 
from  Sorokina  et.  al.  [46],  The  absorption  and  emission  cross-sections  at  the  pump  wavelength 
of  1.9  pm,  opa  and  ape  respectively,  are  0.6  x  10'  cm  and  0.3  x  10'  cm"  while  the  emission 
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cross-section  at  the  input  signal  wavelength  of  2.5  pm,  ase,  is  0.9  x  10'  cm".  The  Cr  doping 
level,  na,  is  calculated  from  the  measured  amount  of  a  small  signal  (i.e.  not  close  to  saturating  or 

2d- 

bleaching  the  absorption)  beam  power  transmitted  through  a  Cr  :ZnSe  sample.  A  typical  value 
of  no  is  1  x  10 19  cm"  . 

The  electron  configuration  for  a  chromium  ion  is  [Ar]3d54s!  and  when  doped  into  a  II -VI 
material  like  ZnSe,  it  will  replace  the  group  II  element,  zinc,  as  Cr"  [55-57].  The  electrons  in 
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the  partially  filled  d  orbitals  interact  strongly  with  the  host  crystal  lattice  giving  rise  to  the  broad 
tunability  of  Cr2+.  The  Cr2+  ion  in  a  tetrahedral  crystal  field  is  a  4-level  system  with  fast 
relaxation  in  the  vibronic  sublevels  which  means  that  nQ  is  effectively  the  sum  of  n*  and  ng, 
where  n*  is  the  upper  state  population  and  ng  is  the  unexcited  (ground  state)  population.  The 
energy  level  diagram  shown  in  Figure  27  shows  the  Cr2+:ZnSe  levels  relative  to  the  band  gap  and 
where  the  ground  state  is  5T2  (ng)  and  the  excited  state  is  5E  (n  ). 


Figure  27:  Cr2+  energy  levels  in  Cr2+:ZnSe  relative  to  the  band  gap  |56] 

In  the  case  of  a  collinearly-pumped  amplifier  where  both  the  pump  and  signal  beams  propagate 
in  the  z-direction,  the  intensity-driven  n  will  decrease  with  z  due  to  the  same  absorption  which 
allows  amplification.  Thus  the  gain  at  each  z  location  will  also  vary  due  to  the  intensity 
dependence  of  the  pumping  rate.  In  order  to  calculate  gain,  G,  in  our  amplifier  we  need  to 
understand  the  relationships  between  n  *,  the  pump  intensity,  Ip  and  the  input  signal  intensity,  Is 
through  a  system  of  differential  equations. 

The  differential  equation  for  n*  will  have  a  positive  term  from  absorption  at  the  pump 
wavelength  and  negative  terms  from  spontaneous  emission  and  stimulated  emission  at  the  signal 
wavelength.  In  addition  to  these,  because  the  absorption  and  emission  bands  of  Cr2  :ZnSe 
overlap  in  the  1.9  pm  region,  there  will  also  be  a  stimulated  re-emission  tenn  at  the  pump 
wavelength  as  well.  The  absorption  cross  section  at  2.5  pm  is  negligible  in  Cr2  :ZnSe  so  there 
will  be  no  absorption  tenn  at  the  signal  wavelength.  The  differential  equation  for  the  pump 
radiation  intensity  will  have  negative  terms  from  absorption  and  passive  losses  at  the  pump 
wavelength  in  addition  to  a  positive  tenn  due  to  stimulated  re-emission  at  the  pump  wavelength. 
The  differential  equation  for  input  signal  intensity  has  a  gain  contribution  from  stimulated 
emission  and  passive  losses  at  the  signal  wavelength. 
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This  system  can  be  described  as 
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where  vp  and  vs  are  the  pump  and  signal  frequencies  and  x  is  the  lifetime  of  the  excited  state.  At 
equilibrium  and  with  no  other  time-dependence,  we  can  solve  (5. 1)  to  get 
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For  a  Cr~  doping  level  of  no,  at  room  temperature  x  is  typically  about  ~4  ps  [58].  It  is  interesting 
to  note  that,  for  a  pump  wavelength  of  1.9  pm, 
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which  in  this  case  means  that  n  can  reach  a  maximum  of  213 no. 

It  is  convenient  to  define  the  pump  and  signal  intensities  in  dimensionless  units  so  that 


r  (  \  7p(Z’r)[  ,  Jp(Z’r) 

JP(Z’r)  =  -h^K«+Vpek  = 


I 


p,sat 


Js(z,r)  = 


JAz’r)„ 


hv,. 


-cr  r  = 


(5.6) 

(5.7) 


where  ISiSat  and  Is,sat  are  the  saturation  intensities  for  signal  and  pump  wavelengths.  With  these 
definitions,  neglecting  radial  dependence  and  assuming  negligible  passive  losses  gives  the 
system 
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These  equations  can  be  solved  numerically  for  some  initial  pump  and  signal  intensities,  Ipo  and 
Iso  and  a  crystal  length  L.  This  system  assumes  that  both  the  signal  and  pump  beams  of  a  given 
power  have  a  uniform  transverse  intensity  over  a  circular  area  A  (with  say  a  radius  of  r).  In 
practice,  however,  the  pump  and  signal  beams  both  have  some  non-unifonn,  mostly  Gaussian 
transverse  intensity  distribution.  In  order  to  model  the  radial  dependence,  we  must  run  the  z- 
propagation  model  for  points  along  the  radius  with  intensities  given  by  the  Gaussian  distribution 
(5. 1 1)  of  r  (with  a  1/e"  radius  of  3  and  a  Gaussian  order  N)  and  integrate  them  over  r  and  9  to  get 
an  output  power. 

l(r,S,N)  =  Qxp  (5.11) 

This  can  take  into  account  non-uniform  saturation  effects  in  the  center  of  the  signal  beam  where 
intensities  are  higher.  The  end  result  of  which  could  be  beams  which  become  less  Gaussian  as 
they  are  amplified. 

A  working  model  was  developed  in  MATLAB  to  represent  our  desired  experimental 
configuration.  Example  model  predictions  of  amplifier  perfonnance  as  a  function  of  varying 
input  signal  levels  and  spot  sizes  can  be  seen  Figure  28,  below,  for  a  crystal  length  of  7.8  mm. 
The  gain  shows  little  sign  of  saturation  effects  due  to  pump  or  signal  intensities  until  very  high 
values  (~  1 0 1 0  W/m2).  The  information  gained  from  the  model  results  was  used  in  the  next 
section  to  drive  the  experimental  design  of  the  MOPA  system.  From  these  models,  the  amplifier 
efficiency  for  a  typical  beam  size  of  40-80  pm  radius  is  around  50%,  which  is  similar  to  the 
types  of  efficiencies  we  see  in  high  power  laser  oscillators.  Therefore,  as  long  as  good  extraction 
efficiency  can  be  achieved  experimentally,  this  approach  is  very  promising  for  power  scaling  of 
Cr24  :ZnSe  laser  systems. 
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Figure  28:  Example  model  predictions  showing  (left)  input  signal  growth  and  pump  absorption  over  the 
crystal  length  for  5  W  input  signal  and  5  W  pump,  both  at  80  pm  spot  size  and  (right)  amplification 
comparisons  for  5  and  10  W  input  signal  powers  with  varying  spot  size  as  a  function  of  pump  power. 

5.3  Experimental  Design 

One  of  the  advantages  of  the  MOPA  design  is  the  flexibility  to  use  almost  any  combination  of 
oscillator  and  amplifier.  A  variety  of  MO  designs  were  built  and  tested  with  the  simplest  PA 
design.  This  was  done  to  evaluate  the  advantages  of  each  design  while  adhering  to  the  optimum 
PA  design  with  respect  to  simplicity  and  perfonnance  predicted  by  the  model.  MO  designs  were 
selected  for  their  ability  to  optimize  certain  operating  modes  such  as  tunability,  narrow  linewidth 
or  power  output. 

There  are  a  number  of  situations  which  may  require  the  orientation  of  the  laser  crystal  at 
Brewster’s  angle  in  a  laser  cavity.  Anti-reflective  (AR)  coatings  in  the  2-3  pm  region  may  be 
too  costly  or  not  robust  enough  for  use  and  thus  another  method  is  needed  to  reduce  reflective 
losses.  However,  using  a  crystal  at  Brewster’s  angle  carries  with  it  some  negative  astigmatism 
effects  which  are  compounded  by  thermal  lensing  in  the  crystal  which  is  now  at  an  angle.  To 
counteract  this,  advantage  can  be  taken  of  the  astigmatism-compensating  properties  of  cavity 
configurations  first  used  by  Kogelnik  in  dye  laser  cavities  [59]. 
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Figure  29:  Astigmatically-compensated  Z-cavity  design  which  can  use  either  a  single  collimated-leg  design 
(dashed)  or  a  double-collimated  leg  design.  OC  =  output  coupler,  LP  =  laser  pump  lens,  FE  =  functional 
element. 

Our  initial  approach  used  a  Kogelnik-type  cavity  with  a  single  collimated  leg  represented  in 
Figure  29  by  the  dashed  lines.  The  cavity  consisted  of  a  flat  45°  dichroic  input  coupler  (M2), 
AR-coated  for  1.9  pm  and  high  reflectivity  (FIR)-coated  for  2-3  pm;  two  5  cm  radius  of 
curvature  (ROC)  spherical  mirrors  (Ml);  and  a  flat  CaF2  50%R  output  coupler  (OC)  located  ~16 
cm  from  Ml.  The  pump  focusing  lens  (LP)  had  an  experimentally-determined  optimal  focal 
length  of  5  cm.  The  polycrystalline,  8. 9x7. 0x3.0  mm  (LxWxH)  Cr  :ZnSe  gain  element  was 
pumped  along  the  long  dimension  with  a  Tm-fiber  laser  linearly  polarized  parallel  to  the  plane  of 
the  diagram.  The  fiber  nature  of  the  pump  system  and  the  lack  of  high-power  isolators  at  1.9  pm 
required  that  this  and  all  subsequent  designs  incorporate  a  non-nonnal  incidence  input  coupler  to 
prevent  back-reflection.  This  cavity  was  shown  to  be  capable  of  over  12  W  of  pure  CW  power 
[4]  and  better  than  40%  slope  efficiency.  This  cavity  was  limited  however,  in  that  its  output 
beam  quality  and  stability  were  poor  at  higher  operating  powers  which  essentially  precluded 
quality  focusing  and  overlap  in  an  amplifier.  It  also  had  no  good  location  in  which  to  insert 
functional  elements  such  as  a  grating  for  tunability  or  a  saturable  absorber  for  modelocking.  As 
can  be  seen  in  Figure  30,  this  cavity  has  periods  of  instability  in  its  output  (denoted  by  larger 
error  bars).  This  was  postulated  to  be  due  to  the  cavity  changing  modes  to  compensate  for 
increasing  thermal  load  in  the  system. 
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Figure  30:  Output  power  vs  incident  pump  power  for  single-leg  Z-cavity.  The  average  slope  efficiency  is 
-40%.  The  error  bars  show  the  minimum  and  maximum  of  1000  measurements  of  output  power  per  second 
at  each  data  point.  The  output  power  oscillation  and  random  non-linearity  of  the  curve  are  attributed  to 
thermal  lensing  instability  and  the  random  pump  beam  distortions  caused  by  atmospheric  water  vapor 
absorption  of  the  pump  radiation. 

To  address  these  concerns,  a  double  collimated-leg  Z-cavity  was  assembled  as  illustrated  with 
the  solid  elements  in  Figure  29.  The  same  5  cm  ROC  mirror  (Ml)  used  as  an  end  mirror  in  the 
above  cavity  was  moved  one  focal  length  from  the  gain  element  and  used  as  the  input  coupler 
while  also  creating  the  second  collimated  leg.  The  ~16  cm  long  ‘functional  element’  (FE)  leg 
was  terminated  with  a  planar  HR  mirror  for  free-running  mode  or  a  300  gr/mm,  2.5  pm  blazed 
grating  in  the  Littrow  configuration  for  tuning  and  line  narrowing.  It  was  found  experimentally 
that  the  optimal  pump  lens  (LP)  focal  length  was  7.5  cm.  Modeling  cavity  mode  overlap  with 
the  pump  can  be  accomplished  using  a  standard  ABCD  matrix  method,  but  the  results  were 
merely  a  good  starting  point  experimentally  due  to  the  variable  and  strong  nature  of  the  thermal 
lensing. 

When  robust  AR  coatings  are  available,  some  of  the  astigmatism  problems  can  be  eliminated  by 
utilizing  a  normal-incidence  Z-cavity  pumping  scheme.  This  simplifies  the  cavity  considerably 
as  seen  in  Figure  31  A.  Cavity  stability  can  be  achieved  using  either  a  concave  HR  end  mirror 
(Za-cavity)  or  a  planar  HR  end  mirror  with  an  intracavity  AR  coated  positive  lens  (Z/,-cavity). 
The  latter  configuration  is  more  complicated  to  construct  in  confined  quarters  and  may  have 
some  loss  depending  on  the  quality  of  the  AR  coating  but  M2  could  be  replaced  with  a  grating  in 
the  Littrow  configuration  to  achieve  cavity  tuning. 
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Figure  31:  A)  Normal-incidence  Z-cavity  scheme  which  can  use  either  an  HR  curved  mirror  (Ml)  or  a 
combination  of  planar  HR  (M2)  and  positive  lens  (LI).  B)  Power  amplifier  design  where  LMO  and  LPA  are 
chosen  for  optimal  overlap  of  pump  and  signal  beams. 

The  power  amplifier  section  used  a  co-linearly  pumped,  polycrystalline,  8. 9x7. 0x3.0  mm 
(LxWxH)  Cr2+:ZnSe  crystal  similar  to  the  MO  crystal,  but  AR  coated  for  1.9-3  pm.  As  seen  in 
Figure  3  IB,  the  150  W  maximum,  1.9  pm,  Tm- fiber  laser  pump  beam  is  co-aligned  with  the  MO 
output.  Both  beams  are  focused  for  optimal  mode  overlap  of  around  80  pm  radius  with  15  cm 
focal  length  lenses  (LMo  and  LPa). 

5.4  Experimental  Results 

The  Z-  and  Z-cavity  oscillators  were  optimized  for  highest  output  at  the  upper  limit  of  their 
stability  region  (Z-cavity)  or  just  below  the  damage  threshold  for  the  AR  coatings  (Z-cavities). 
All  experiments  were  done  with  the  laser  samples  cooled  with  circulating  water  maintained  at 
10°  C.  Experimentally,  it  was  determined  that  the  Lb- cavity  was  easily  out-perfonned  by  the  La- 
cavity,  possibly  due  to  additional  losses  from  the  intra-cavity  lens.  Once  an  MO  was  optimized 
in  this  way,  the  PA  pump  power  was  then  varied  to  test  the  system  and  compare  to  the  theoretical 
model.  The  data  shown  in  Figure  32  represents  the  best  performance  of  each  of  the  MOPA 
configurations.  From  this,  it  can  be  seen  that  the  La- cavity  is  better  suited  for  power  scaling  than 
the  Z-cavity.  This  is  likely  due  to  the  limited  ability  of  the  fixed-angle  curved  mirrors  in  the  Z- 
cavity  to  compensate  for  the  increasing  astigmatism  from  the  Brewster-angle  thermal  lens  in  the 
gain  element.  However,  it  should  be  noted  that  the  output  from  the  Z-cavity  is  tunable  and,  due 
to  the  poor  performance  of  the  Z/, -cavity,  this  is  not  an  option  in  the  linear  cavity  configuration. 
The  experimental  data  follows  the  model-predicted  performance  at  first  but  then  under-performs 
the  predictions.  This  could  be  due  to  thermal  lensing  in  the  amplifier  creating  non-optimal 
overlap  between  the  pump  and  signal.  Both  signal  and  pump  should,  in  general,  experience  the 
same  thennal  lensing  but  dispersion  could  result  in  small  differences  in  lensing  power. 
Experimentally,  it  was  noticed  that  any  small  deviations  in  overlap  resulted  in  significant  output 
reduction,  so  even  minute  lensing  variations  could  result  in  the  observed  differences  between 
model  and  experiment. 
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Figure  32:  Power  scaling  ability  of  MOP  A  configurations.  Experimental  data  (discrete  points)  is  compared 
with  model  predictions  for  La-  and  Z-cavities,  both  with  100  pm  spot  sizes. 

The  output  spectrum  of  the  free-running  (no  grating  present)  Lb- cavity  MO  was  measured  as  a 
function  of  output  power.  Spectra  were  measured  by  a  2  MHz  Cal  Sensors  256  element  PbSe 
linear  detector  array  coupled  to  the  exit  slit  of  a  750  mm  Acton  Research  Corporation 
monochromator.  This  system  was  calibrated  and  optimized  using  a  Yokogawa  735305  optical 
spectrum  analyzer  (OS A)  and  a  150  W  IPG  Tm  fiber  laser.  The  OSA  showed  a  0.2  mn  linewidth 
for  the  150  W  IPG  laser  and  after  proper  positioning  of  the  array  at  the  focal  plane  and 
adjustment  of  the  monochromator  slit  widths,  similar  resolution  was  achieved  with  the 
spectrometer  system.  Figure  33  shows  that  the  width  of  the  collection  of  output  peaks  of  the 
laser  is  not  a  function  of  pump  power;  however,  the  weighted  average  wavelength  increases  with 
pump  power,  going  from  2415  to  2442  mn.  There  are  a  number  of  thermal  effects  which  could 
be  postulated  to  shift  the  average  output  wavelength  in  either  direction,  but  it  is  difficult  to 
predict  which  will  dominate.  A  shift  to  the  longer  wavelengths  could  be  due  to  higher 
temperatures  in  the  crystal  affecting  non-radiative  transitions  for  higher  energy  levels. 
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Figure  33:  Output  spectrum  of  free-running  Tft-cavity  as  a  function  of  output  power  showing  a  shift  to  higher 
wavelengths  at  increased  power. 

Representative  tunable  output  characteristics  of  the  amplified,  grating  equipped  Z-cavity  are 
shown  in  Fig.  34.  A  maximum  power  of  8.9  W  was  achieved  at  2450  nm  with  tunable  lasing 
over  the  2275-2700  nm  range,  which  was  limited  by  grating  efficiency  and  wavelength- 
dependent  coating  reflectivity.  Figure  34  shows  the  dependence  of  the  measured  output  power 
on  wavelength  superimposed  with  the  linewidths  measured  at  the  same  points,  which  were  on  the 
order  of  2  nm. 


Figure  34:  Tunable  output  power  of  Z-cavity  with  representative  linewidth  of  output. 


One  of  the  goals  of  the  MOPA  configuration  is  an  improvement  in  the  beam  quality  for  a  given 
output  power.  The  need  for  this  can  be  understood  by  looking  at  the  near-field  beam  profile  of 
the  MO  output.  Figure  35  shows  a  single-frame  excerpt  from  a  video  showing  the  near-field 
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beam  profile  of  the  Z^-cavity  as  a  function  of  decreasing  pump  power.  The  cavity  will  run  stable 
in  one  mode  (TEM0o)  at  low  powers  until  the  power  changes  enough  to  allow  or  cause  a  run  in 
combinations  of  modes  with  the  best  stability.  This  type  of  transition  period  in  the  cavity  is 
marked  by  unstable  and  fluctuating  output  powers  and  beam  quality.  On  the  other  hand,  with  a 
MOPA  system,  the  MO  could  be  placed  in  a  stable  region  while  the  system  power  output  is 
varied  by  the  PA,  allowing  for  more  stable  overall  output  power  and  quality. 


Pump  pow  er  (W) 


Figure  35:  Single-frame  excerpt  (left)  from  video  showing  near-field  beam  profile  changes  of  Lb- cavity  as  a 
function  of  decreasing  pump  power.  Movie  starts  with  maximum  pump  power  and  ends  at  zero  pump  power. 
Transitions  between  different  modes  cause  instability  in  power  output  (right). 


The  ability  of  the  MOPA  system  to  achieve  higher  output  powers  while  maintaining  beam 
quality  was  tested  by  measuring  the  M“  of  the  La-cavity.  At  1  W,  the  cavity  had  an  M  of  1 .2  and 
at  3  W  it  had  an  M2  of  1.8.  With  the  MO  running  at  1  W  and  amplified  to  3  W  with  the  PA,  the 
system  had  an  M"  of  1.3.  This  supports  the  claim  that  a  MOPA  system  can  be  used  to  power 
scale  a  system  where  other  issues  may  limit  oscillator  operation. 

5.5  Conclusions 

2”b 

In  conclusion,  we  have  demonstrated  a  CW  Cr  :ZnSc  MOPA  laser  with  output  power  over  10 
W  and  amplifier  gain  of  2X.  We  developed  an  amplifier  model  for  co-linearly  pumped 
Cr“  :ZnSe,  which  compares  well  with  our  experimental  data.  Use  of  a  MOPA  configuration 
allowed  us  to  demonstrate  high-power  and  good  beam  quality  in  a  Cr2+:ZnSe  laser.  In  addition, 
we  demonstrated  single-knob  tunability  at  multi-watt  powers  over  a  400  nm  tuning  range. 


As  is  common  with  chromium-doped  II- VI  lasers,  thermal  effects  limited  the  output  power  and 
beam  quality  of  this  setup.  Improvements  in  cooling  are  needed  to  take  full  advantage  of  the 
available  pump  power  and  increase  system  stability.  This  MOPA  system  will  be  evaluated  with 
addition  MO  configurations  to  study  the  optimal  operating  conditions  and  increase  output  laser 
power.  Thermal  management  techniques  will  continue  to  be  improved  in  an  effort  to  reduce 
thermal  instability  in  the  oscillator. 
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Chapter  6 
Conclusions 

6.1  Summary 

In  this  research  a  number  of  approaches  were  used  to  understand  and  mitigate  thennal  issues 
limiting  the  level  to  which  tunable  Cr2+:ZnSe  lasers  could  be  power  scaled.  It  was  found  that 
careful  consideration  of  heat  removal,  lensing  compensation  and  system  architecture  can  help 
remove  many  of  the  hurdles  to  high-power  operation.  A  finite-element  analysis  model  was  used 
to  help  understand  the  time-evolution  of  the  thermal  lensing  problem.  Crucial  initial  work  in  the 
development  of  fiber-based  Cr“  :ZnSe  lasers  was  completed  and,  while  not  successful  in 
producing  a  Cr2  :ZnSe  fiber,  helped  understand  the  material  problems  associated  with  this 
geometry.  A  Cr  :ZnSe  laser  was  constructed  which  demonstrated  over  12  W  of  pure  CW  output 
power,  doubling  the  previous  record.  A  model  for  amplification  of  co-linearly  pumped  master 
oscillator  power  amplifier  systems  was  developed  and  used  to  design  the  first  demonstrated 
Cr2+:ZnSe  MOPA  systems  with  greater  than  10  W  pure  CW  power  and  9  W  tunable  over  400 
mn. 

6.2  Future  Work 

This  work  showed  that  there  are  still  a  number  of  things  that  can  be  done  to  power  scale  these 
lasers.  Using  thennal- gradient  reducing  geometries  like  zig-zag  slab  lasers  is  planned  to  further 
increase  output  powers.  The  fiber  laser  work  has  lead  to  a  better  understanding  of  the  difficulties 
involved  in  producing  a  Cr“  :ZnSe  fiber  laser.  This  has  lead  into  new  efforts  to  produce  these 
types  of  waveguides  which  hold  promise.  The  MOPA  system  design  will  also  be  used  to 
produce  high-power  modelocked  pulses  for  use  in  pumping  nonlinear  devices. 
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Appendix  1 

Listing  of  MATLAB  Codes 

Listing  1:  thermal_lens.m 

function  f  =  thermal_lens ( fem,  nx,  ny,  nz,  n  window); 

~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k 

%  FUNCTION  thermal  lens (fem,  nx,  [ny] ,  [nz],  n  window) 

%  function  to  do  line  integral  along  z  of  temperature  and  plot  2D  result 
%  fem:  fem  structure  exported  from  FemLab 

%  nx:  number  of  points  in  x  data  array,  usually  51  or  101 

%  ny:  optional,  number  of  points  in  y  array,  defaults  to  nx 

%  nz:  optional,  number  of  points  in  z  array,  defaults  to  nx 

%  n  window:  window  size  in  number  of  points 

~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k  ~k 

if  nargin  <  2 
nx  =  51; 

end 

if  nargin  <  3 
ny  =  nx; 

end 

if  nargin  <  4 
nz  =  nx; 

end 

if  nargin  <  5 

n  window  =  0; 

end 


dndt 

xmax 

=  7  0  e  -  6  ; 

=  postmax (fem. 

'x'); 

xmin 

=  postmin(fem. 

'x'); 

ymax 

=  postmax (fem. 

'  y ' )  ; 

ymin 

=  postmin(fem. 

'  y ' )  ; 

zmax 

=  postmax (fem. 

'  z  '  ) 

zmin 

=  postmin(fem. 

'  z  '  ) 

%  YOU  MUST  SPECIFY  THE  Z  DIM  FOR  MODELS  WITH  WINDOWS  AND  RING  HEATSINKS 
%zmax  =  0.002; 

%zmin  =  -0.001; 

%  find  out  the  baseline  temp  from  the  constants 
TO  =  postmax (fem,  'TO'); 

%  define  dimension  vectors 
xd  =  linspace (xmin,  xmax,  nx)  ; 
yd  =  linspace (ymin,  ymax,  ny)  ; 
zd  =  linspace ( zmin,  zmax,  nz); 

%  make  vectors  into  a  3D  mesh 
[x, y, z ] =meshgrid (xd,  yd,  zd)  ; 

%  get  the  temperature  data  from  the  'fem'  structure  for  the  mesh 
Txyz  =  postinterp ( fem,  ' T '  ,  [x ( : )  '  ;  y ( : )  '  ;  z ( : )  '  ]  )  ; 

%  replace  any  NaNs  with  293 
i  =  find (isnan (Txyz) ) ; 

Txyz (i)  =  TO; 

%  break  the  long  array  into  a  3D  array 
Txyz=reshape (Txyz, size (x) ) ; 

%  find  out  the  length  of  each  section  of  the  array 
[dx,  dy,  dz]  =  size (Txyz); 

%  average  the  temp  in  the  z  dir  and  subtract  baseline  temp  to  get 
%  ave  temp  rise 
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T_z_ave  =  (sum (Txyz, 3) /dz)  -  TO; 

%  surface  plot  the  average  temp  rise 
xdl  =  xd*1000; 
ydl  =  yd* 1000; 
figure ( 1 ) 

surf (xdl,  ydl,  T_z_ave) ; 

set (gcf ,  'Color',  'white'); 

h  =  findobj (gcf,  'Type',  'surface'); 

set (h,  'LineStyle',  ':'); 

axis  tight; 

zlabel ( ' \DeltaT  Ave ' ) ; 

T  x  =  T  z  ave (:, floor (dx/2 ) +1 ) ;  %  pulls  temperature  profile  for  y=0  in  the  x 
direction 

T  x  =  T  x';  %  transposes  the  matrix  to  get  it  in  the  right  orientation 
max  ave  temporise  =  max(T_x) 

%  Focal  Length  Calculation 

%  first  decide  what  data  window  we  are  fitting  to 
if  n^window  >  0 

%  if  a  window  is  given,  create  the  smaller  data  sets 
index  start  =  (nx  -  l)/2  -  n  window  +1; 

index  end  =  (nx  -  l)/2  +  n  window  +  1; 

x  win  =  xd ( index^start : index  end); 

T  x  win  =  T  x (index  start: index  end); 

else 

%  otherwise,  take  just  the  middle  section  of  the  data  and  x  axis  for 
%  20% 

index  start  =  (nx  -  l)/2  -  round (nx*. 2)  +1; 

index  end  =  (nx  -  l)/2  +  round (nx*. 2)  +  1; 

x  win  =  xd(index  startrindex  end); 

T  x  win  =  T  x (index  startrindex  end); 

end 

%  now  fit  a  polynomial  to  the  half-max 
poly_constants  =  polyf it (x  win,T  x  win, 2); 

%  take  the  first  element,  the  xA2  constant 
b  =  poly_constants ( 1 ) ; 

TOVal  =  poly_constants  ( 3 ) ; 

%  plot  the  full  data 
figure (2 ) 
plot (xd, T_x) ; 

%  plot  the  window  data 
hold  on 

plot (x_win, T  x  win, ' k- ' ) ; 

plot (x  win, TOVal  +  b.*x  win . A2 , ' r : ' , ' LineWidth ' , 2 ) ; 
set  (gcf ,  'Color',  'white'); 
axis  tight; 

ylabel ( ' \DeltaT_Ave ' ) ; 
xlabel ( ' Position ' ) ; 
hold  off 

%  calculate  the  lens  power  in  diopters 
f  =  (4.*(zmax  -  zmin) . *dndt . *b) 
disp ( ' diopters ' ) ; 

fl  =  ['lens  power  =  ',  num2str(f),  '  mA-l']; 
text (0, max  ave  temp  rise/2, fl); 

laser  diameter  fit  =  xd (index  end)  -  xd (index  start) 
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Listing  2:  HSIC.m 

function  IC  =  HSIC(xl,  x2,  yl,  y2 ,  thick,  sigma,  order) 

^•k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k 

%  IC  =  HSIC(xl,  x2,  yl,  y2 ,  thick,  sigma,  order) 

%  Calculates  the  integration  constant  needed  to 
%  normalize  the  absorbed  power  density  to  1 

^•k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k  -k 

%  set  up  the  x  and  y  limits 
x  =  linspace(xl,  x2,  101); 

y  =  linspace(yl,  y2,  101); 

%  create  a  2D  mesh  from  those  co-ords 
[x,  y]  =  meshgrid(x,  y)  ; 

%  calculate  the  heat  source  for  each  co-ord  in  the  mesh 
H  =  HSG(x,  y,  thick,  sigma,  order); 

%  integrate  the  2D  array  to  get  the  Absorbed  Power  Density  (volume) 

%  do  not  forget  to  multiply  by  the  thickness  to  account  for 
%  the  integration  in  the  z  direction  which  is  considered 
%  to  be  uniform  for  this  model 

%  note  that  integrations  in  mat-lab  need  to  be  done  using  the  equivalent 
%  of  a  Riemann  Sum:  sum  of  f (x) *dx 
dx  =  (x2  -  xl)/100; 
dy  =  (y2  -  yl) /100; 

APD  =  sum ( sum (H*dx) *dy) *thick; 

IC  =  1 /APD; 

comp  =  2/ (pi*sigmaA2*thick) ; 
ratio  =  IC/comp; 

disp (' integration  constant  =  '); 
disp (ratio) ; 
disp  (  ' *  '  )  ; 

disp ( ' 2/ (pi*sigmaA2*zmax)  ' )  ; 
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Listing  3:  HSG.m 

function  f  =  HSG (x, y, zmax, sigma,  order) 

^•k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%  f  =  HSG (X, Y, ZMAX, SIGMA) 

%  HEAT  SOURCE  GAUSSIAN 
%  x  =  x  co-ord 

%  y  =  y  co-ord 

%  zmax  =  z  thickness 
%  sigma  =  pump  beam  radius 

^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

if  nargin  <  5 
order  =  1; 

end 

xO  =  0 ; 
yO  =  0; 

%  normalized 

%f  =  (2/ (pi*sigma*sigma*zmax) ) *exp (- (2* ( (x  -  xO).A2  +(y- 
yO )  . A2 ) )  / (sigma. A 2) ) ; 

%  unnormalized 

%  any  order  other  than  1  needs  to  be  divided  by  2 
if  order  >  1 

order  =  order/2; 

end 

f  =  exp(-(2*((x  -  xO).A2  +  (y  -  yO ). A2 ). "order) /( sigma . A (2*order) ) ) 
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Listing  4:  HSICR.m 

function  IC  =  HSICR(rmax,  thick,  sigma,  order) 

-k  -k  -k  -k  -k  -k  -k  -k  -k  -k 

%  IC  =  HSICR(r,  thick,  sigma,  order) 

%  Calculates  the  integration  constant  needed  to 
%  normalize  the  absorbed  power  density  to  1 
%  for  a  radial  system  =  round  disk 

%rk  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 
o 

%  set  up  the  x  and  y  limits 
r  =  linspace(0,  rmax,  1001); 

%  calculate  the  heat  source  for  each  co-ord  along  the  radius 
H  =  HSGR(r,  sigma,  order); 

%  integrate  to  get  the  Absorbed  Power  Density  (volume) 

%  do  not  forget  to  multiply  by  the  thickness  to  account  for 
%  the  integration  in  the  z  direction  which  is  considered 
%  to  be  uniform  for  this  model 
%  and  also  by  2*pi  for  the  radial  symmetry 

%  note  that  integrations  in  mat-lab  need  to  be  done  using  the  equivalent 
%  of  a  Riemann  Sum:  sum  of  f (x) *dx 
dr  =  rmax/1000; 

APD  =  sum (H . *r*dr ) * (2*pi ) *thick; 

IC  =  1/APD; 

comp  =  2/ (pi*sigmaA2*thick) ; 
ratio  =  IC/comp; 

disp (' integration  constant  =  '); 
disp (ratio) ; 
disp  (  ' *  '  )  ; 

disp ( ' 2/ (pi*sigmaA2*zmax)  ' )  ; 
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Listing  5:  HSGR.m 

function  f  =  HSGR(r,  sigma,  order,  rO) 

S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

%  f  =  HSGR (r, SIGMA,  ORDER) 

%  HEAT  SOURCE  GAUSSIAN  -  radial 
%  r  =  radial  co-ord 
%  sigma  =  pump  beam  radius 

%  order  =  gaussian  order 

S^k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k  k 

if  nargin  <  3 
order  =  1; 

end 

if  nargin  <  4 
rO  =  0; 

end 

%  unnormalized 

%  any  order  other  than  1  needs  to  be  divided  by  2 
if  order  >  1 

order  =  order/2; 

end 

f  =  exp ( - (2* ( (r  -  rO).A2). ''order) / (sigma . A (2*order ) ) ) 
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Listing  5:  gain mo del3. m 

function  g  of  Pp  =  gain  model3(Pp,  PsO,r,L  s,n  0,tau,nz,nr) 

%  function  g_of_Pp  =  gain  model2(Pp,  PsO,  r,  L_s, n_0, tau, nz, nr) 

np  =  size (Pp, 2) ; 
g_of_Pp  =  zeros (l,np); 

%  create  r  variables 
if  nargin  <  8 
nr  =  101; 

end 

if  nargin  <  7 
n  z  =  101; 

end 

if  nargin  <  6 

tau  =  4e-6;  %lifetime  in  s 

end 

if  nargin  <  5 

n  0  =  10e24;  %  l/mA3  doping  level 

end 

if  nargin  <  4 

L  s  =  7.8e-3;  %m  sample  length 

end 

if  nargin  <  3 

r  =  80e-6;  %m  spot  radius 

end 

rp  =  linspace ( 0 , 3*r , nr ) ; 
rs  =  rp; 

%  define  intensities  I  (r) 

Irp  =  HSGR(rp,  r,  1); 

Irs  =  Irp; 

drp  =  rp  ( 2  )  -  rp  ( 1 )  ; 

drs  =  rs(2)-rs(l); 

%  pre-allocate  vars  for  speed 
Ps_out  =  zeros (np, 1) ; 

Pp_out  =  zeros (np, 1) ; 

Is_out  =  zeros (1, nr); 

Ip_out  =  zeros (1, nr); 

gj  =  zeros (1, nr);  %  gain(j) 

paj  =  zeros (1, nr);  %  pump  absorbed 

%  calculate  non-normalized  power 

Ppn  =  sum ( Irp . *rp*drp) *2*pi ; 

Psn  =  sum ( Irs . *rs*drs ) *2 *pi ; 

%  calculate  int  const 
ICs  =  PsO/Psn; 

%  now  define  normalized  intensity  signal 
Inrs  =  ICs*Irs; 

%f or  i=l : np 

%  [gl,  pi]  =  amplif ier_gain3 (Pp (i) /pi/rA2, PsO/pi/rA2, L_s, n_0 , tau, nz) 
%  Ps_out(i)  =  gl*PsO; 

%  Pp_out(i)  =  pl*Pp(i); 

%end 

%plot (Pp, Ps_out,  Pp,  Pp_out) ; 
for  i=l : np 

%  amp  int  constant 
ICa  =  Pp (i) /Ppn; 

%  normalized  amp  intensity 
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Inrp  =  ICa*Irp; 
for  j  =  1 : nr 

[gj  (1, j ) ,  paj  (1, j  )  ]  = 

amplifier  gain3 ( Inrp ( j ) , Inrs ( j ) , L  s,n  0,tau,nz) 
end 

Is_out  =  gj.*Inrs; 

Ip_out  =  pa j.* Inrp; 

Ps_out(i)  =  surafls  out . *rs*drs) *2*pi; 

Pp  out(i)  =  sum(Ip  out.*rp*drp)*2*pi; 
g_of_Pp(l,i)  =  Ps_out ( i ) /PsO ; 

end 
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Listing  6:  amplifier  gain3.m 

function  [g,  pabs]  =  amplifier  gain3 (IpO, IsO, L_s, n_0, tau,  nz) 
%  this  function  will  calculate  the  gain  of  an  amplfier 
%  using  numerical  model 


lambda  s  =  2.5e-6;  %m  signal  wavelength 
lambda_p  =  1.9e-6;  %m  pump  wavelength 
nu  s  =  cO/lambda  s;  %1/s  signal  freq 
nu_p  =  cO/lambda_p;  %1/s  pump  freq 
h  =  6 . 626068e-34;  %mA2kg/s  planck 


pump  absorption  xsect 
pump  emission  xsect 
signal  emission  xsect 


^lifetime  in  s 


l/mA3  doping  level 


sig_pa  =  6e-23;  %mA2 
sig_pe  =  3e-23;  %mA2 
sigse  =  9e-23;  %mA2 
if  nargin  <  6 
n  z  =  101; 

end 

if  nargin  <  5 
tau  =  4e-6 ; 

end 

if  nargin  <  4 

nj]  =  10e24 ; 

end 

if  nargin  <  3 

L_s  =  7.8e-3; 

end 

%Is_sat  =  h*nu_s/ (sig_pe*tau) ;  %W/mA2 
%Ip_sat  =  h*nu_p/ (sig_pa+sig_pe) /tau; 
z  =  linspace ( 0 , L_s , nz ) ;  %sample  coords 
dz  =  z  (2)  -  z  (1)  ; 

%  pre-allocate  some  variables 
gz  =  zeros ( 1 , nz ) ; 

Jpz  =  zeros (l,nz) ; 

Jsz  =  zeros (l,nz); 

Jpz(l)  =  IpO* (sig_pa+sig_pe) *tau/ (h*nu_p) ; 

Jsz(l)  =  IsO*sig_se*tau/ (h*nu_s) ; 

gz ( 1 )  =  1; 

k_0  =  sig_pa*n_0; 


fem  sample  length 


sig  saturation  inten 
%pump  sat  inten 


for  int  =  2 : nz 

gz(int)  =  k_0* (sig_se/ (sig_pa+sig_pe) ) * ( Jpz (int-1) / ( Jpz (int-1) 
1)  +  1)  )  ; 

Jsz (int)  =  Jsz (int-1)  +  Jsz (int-1) *gz (int) *dz; 

Jpz (int)  =  Jpz (int-1)  +  Jpz (int-1 )* (-k  0  + 

(sig_pa+sig_pe) *gz (int) /sig_se) *dz; 
end 

g  =  Jsz (nz) /Jsz (1) ; 
pabs  =  Jpz (nz ) / Jpz ( 1 )  ; 


+  Jsz (int- 
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Appendix  2 

Lab  View  Control  Program  for  Laser  and  Power  Meter  Control 


A  revised  control  program  was  written  and  tested  for  use  with  both  IPG  lasers  and  all  Newport 
power  meters  (1918-C  and  2935  models).  This  program  allows  for  simultaneous  control  of  both 
lasers  independently  while  collecting  statistical  data  from  the  power  meter.  Data  collected  from 
the  power  meter  can  be  saved  for  use  in  MatLab  or  other  data  processing  program.  The  status  of 
both  lasers  is  displayed  in  real  time  to  alert  the  user  to  any  problems. 

The  LabView  programming  environment  is  an  icon-based  system.  Typically  programs  are 
meant  to  execute  once  or  continuously  in  a  loop.  For  the  command-based  terminal  language  of 
the  both  the  IPG  lasers  and  the  Newport  power  meters,  an  even-driven  architecture  is  more 
suitable  so  that  queries  are  not  constantly  being  sent  to  all  instruments.  To  accomplish  this  with 
LabView,  it  is  necessary  to  use  the  LabView  event  structure  to  respond  to  button  pushes  and 
similar  events. 

The  150  W  laser  (Figure  36)  is  controlled  by  entering  changes  to  percent  current  (from  13  to 
100%).  The  40  W  (Figure  37)  laser  is  controlled  by  entering  either  changes  to  power  level  (from 
1-40  W)  or  changes  to  current  (0 — 5000  inA).  The  status  display  block  for  either  will  show  any 
status  changes  reported  by  the  laser. 


IPG  150  W 

Select  Laser(s) 
Before  Run 


STOP 


150  W  Com  Port 
frlCOMl  - 

Emission  Align  Beam 


Set  Current  %  Enter 


rllO.OO 


~OK~l 


Temperature 


Current  Setpoint  % 


Power  Output 


Firmware 


Refresh  Status 
|  OK  | 


)  Overtemperature 

)  Emission  on  or  starting 

)  High  backreflection!! 

)  Analog  control  enabled 

^  Fiber  fuse  detected 

^  Module(s)  disconnected 

^  Module(s)  failed 

9  Aiming  beam  on 

Interlock  open 

^  Optical  interlock  failed 

^  Power  supply  off 

®  Modulation  enabled 


Figure  36:  Control  panel  for  IPG  150W  Tm:Fiber  laser  control 
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IPG  40  W 


STOP 


40  W  Com  Port 


Temperature  2 


<6|COM2 

1 

Set  Power  (W)  Set  P 

|  ok  1 

Set  Current  (mA)  Set  I 

^jo.oo 


OK~l 


Current  (mA) 


Power  Output  2 


Firmware  2 


)  Power  Mode 
)  High  backreflection!! 

|  Interlock  open 
)  Overheating 

> 

|  Emission  on  or  starting 
Refresh  Status  2 
1  OK  | 


Figure  37:  Control  panel  for  IPG  40W  Tm:Fiber  laser  control 


The  control  block  for  the  power  meter  (Figure  38)  allows  for  control  of  either  a  1918-C  single 
channel  power  meter  or  a  2935  dual-channel  power  meter.  The  control  panel  can  be  used  to  zero 
either  channel  and  take  set  numbers  of  data  points  at  a  given  rate  and  then  store  the  mean,  min, 
max  and  standard  deviation.  It  should  be  noted  that  the  rate  is  in  approximately  kHz  (i.e.  10000 
pts  at  “1”  equals  about  10  seconds).  The  right  side  of  this  panel  can  be  used  to  run  automated 
scans  by  starting  one  of  the  lasers  at  a  given  current  (again  %  for  the  150  and  mA  for  the  40), 
waiting  a  set  amount  of  time,  then  taking  a  given  number  data  points,  stepping  by  a  given  step, 
until  reaching  the  ending  current.  This  is  useful  for  automated  slope  efficiency  measurements 
which  remove  human  variations  from  the  results. 
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Power  Meter  Channel  Zero  Detector  Model 

Brightness  B  powe|-  |  Run  | 

$  |o 


Mean  Power 

Max. 

& 

v — ml 

Std.  Dev. 

Minimum 

r 

15 

Collect  Stats 

1  OK  | 

Delete  Row  Save  Data  Collect  Single  Chan 


Start  Collection 
1  OK  | 


Start  (mA) 


I  c 


0 


40W 


Figure  38:  Control  panel  for  power  meter  automated  data  collection  (left)  and  for  automated  laser  control 
with  data  collection  (right) 
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Appendix  3 

Z-Cavity  Alignment  Procedure 

The  procedure  described  below  is  one  used  to  align  the  single-leg  Z-cavity  described  in  Chapter 
5.  This  same  procedure  was  adapted  for  use  with  the  other  MO  cavities.  In  Figure  39  below, 
Ml,  M2  are  50  mm  ROC  mirrors,  DC  is  a  flat  in-coupler,  all  of  which  are  HR@2.5  pm,  AR  @ 
1. 9-2.1.  LI  is  a  50  mm  focal  length  lens,  dl,  d2  and  d3  are  25mm,  d4  is  160  mm,  and  A1-A5  are 
adjustable  aperture  irises.  The  1.907  pm  pump  comes  into  LI  co-aligned  with  a  5  mW  Helium 
Neon  alignment  laser.  LI,  Ml  and  M2  are  on  translation  stages. 


Figure  39:  High  power  Kogelnik  cavity 


The  following  rigorous  alignment  procedure  was  developed  with  input  from  collaborators  at 
UAB  and  used  to  align  the  laser  cavity: 

1.  Remove  all  optical  elements  from  the  optical  table.  Ensure  that  the  HeNe  beam  is 
parallel  to  the  optical  table  and  is  aligned  with  the  pump  beam. 

2.  Install  the  input  coupler  (DC). 

3.  Position  an  iris  aperture  (Al)  -150  mm  from  the  proposed  location  of  the  pump  lens 
and  center  it  with  the  HeNe  beam. 

4.  Add  another  iris  aperture  (A2)  before  the  pump  lens  position  and  center  it  with  the 
HeNe  beam. 

5.  Install  the  pump  lens  (LI)  and  align  it  vertically  and  horizontally  to  make  sure  that  the 
propagated  HeNe  beam  is  centered  on  the  iris  aperture  Al. 
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6.  Install  an  iris  aperture  (A5)  and  center  it  with  the  reflection  of  the  HeNe  beam  from  the 
internal  surface  of  the  flat  pump  mirror.  The  in-coupler  DC  causes  significant 
diffraction.  The  iris  should  be  aligned  with  the  second  brightest  reflected  spot. 

7.  Install  the  laser  crystal  and  align  it  to  the  Brewster  angle  (67  degrees  for  zinc  selenide). 
This  can  be  done  by  using  a  polarizer  on  the  HeNe  and  finding  the  minimum  reflection 
point.  Triangulation  can  also  be  used.  The  laser  crystal  should  be  mounted  on  a  three 
axis  stage. 

8.  Even  a  very  weak  HeNe  beam  should  be  visible  after  propagation  through  the  crystal, 
so  ensure  that  the  beam  appears  somewhere  on  iris  aperture  Al.  It  will  be  displaced 
from  passing  through  the  crystal. 

9.  Install  mirror  M2  and  adjust  it  to  the  appropriate  angle  (15-20°).  Displace  the  mirror 
vertically  as  necessary  until  the  reflected  HeNe  beam  d4  is  horizontal. 

10.  Adjust  the  distance  d3  until  the  reflected  d4  HeNe  beam  is  collimated.  This  indicates 
optimal  distance  for  d3. 

11.  Use  two  iris  apertures  (A3,  A4)  to  indicate  the  path  of  the  d4  beam  leg.  Care  should 
be  taken  with  these  apertures  to  not  interfere  with  the  cavity  mode  while  leaving  room 
for  the  outcoupler  OC. 

12.  Align  a  second  HeNe  beam  counter-propagating  along  d4  through  apertures  A3  and 
A4.  If  the  reverse  HeNe  beam  is  aligned  correctly  with  d4,  you  should  also  observe  its 
collimated  spot  centered  at  the  second  iris  aperture  A2. 

13.  Install  the  end  mirror  Ml  and  align  it  so  that  the  reversed  HeNe  beam,  reflected  from 
Ml  is  centered  at  the  aperture  A5.  This  is  the  hardest  spot  to  align  in  the  procedure. 
Turning  out  the  room  lights  is  helpful. 

14.  With  no  outcoupler,  set  up  a  suitable  detector  (like  an  extended  range  InGaAs 
detector)  and  a  chopper  at  the  output  end  of  the  propagation  path,  past  A4  and 
off/through  a  filtering  dichroic  mirror  FI  (to  dump  any  pump  signal).  Turn  on  the 
pump  laser  and  set  it  to  2-3  W.  Observe  the  fluorescence  signal  on  an  oscilloscope 
triggered  with  the  chopper  signal. 

15.  Maximize  the  fluorescence  signal  at  the  detector  by  longitudinal  adjustment  of  the 
pump  lens  and  fold  mirror  M2. 

16.  Maximize  the  fluorescence  signal  at  the  detector  by  angular  alignment  of  the  end 
mirror  Ml. 

17.  Turn  off  the  pump. 

18.  Install  the  output  coupler.  Use  a  highly  reflective  OC  (95%)  for  easier  initial 
alignment.  Use  the  retro-reflection  of  the  reverse  HeNe  for  preliminary  alignment. 

19.  There  is  no  step  19. 

20.  If  the  forward  HeNe  beam  reflection  from  OC  is  visible  through  the  cavity,  align  OC 
so  that  the  reflected  beam  passes  through  the  aperture  A3. 

2 1 .  Turn  the  pump  laser  on  and  set  it  to  5  W. 

22.  Quickly  scan  the  OC  horizontally  while  slowly  scanning  it  vertically  until  a  lasing 
spike  is  seen  in  the  fluorescence  signal. 

23.  Further  align  all  other  optical  elements  in  a  repeatable  pattern,  slowly  decreasing  the 
pump  power  after  each  alignment  cycle,  until  the  minimum  lasing  threshold  is  obtained. 

24.  If  lasing  could  not  be  obtained  on  the  previous  step  then  the  cavity  is  probably  in 
unstable  condition. 

25.  Turn  the  pump  laser  off  and  pre-align  the  OC  again  as  in  steps  18-20. 
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26.  Turn  on  the  pump  and  set  it  to  5  W. 

27.  Move  the  fold  mirror  M2  closer  to  the  crystal  by  a  few  millimeters. 

28.  Slowly  increase  d3  in  small  steps  (-200-500  pm)  and  at  each  step  shake  the  OC  in 
two  angular  directions.  The  stable  alignment  of  the  laser  cavity  will  eventually  be  found 
and  lasing  obtained. 

29.  Optimize  the  cavity  to  the  minimum  threshold  by  realigning  all  optical  elements  in  a 
repeatable  pattern. 

30.  Note  that  the  distance  d4  is  not  critical.  It  can  vary  from  30  mm  to  100  cm  without 
altering  the  cavity  stability. 
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